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Abstract 

Addition tlieorems are principal tools that express a function f{r ± r') in terms of products of other functions 
that only depend on either r or r'. The best known example of such an addition theorem is the Laplace expansion 
of the Coulomb potential which possesses a characteristic two-range form. Guseinov [Chem. Phys. 309, 209 - 213 
(2005)] derived one-range addition theorems for the Coulomb potential via the limit /3 ^ in previously derived 
one-range addition theorems for the Yukawa potential exp(^—(3\r — r'l)/|r — r'\. At first sight, this looks like a 
remarkable achievement, but from a mathematical point of view, Guseinov's work is at best questionable and in some 
cases fundamentally flawed. One-range addition theorems are expansions in terms of functions that are complete and 
orthonormal in a given Hilbert space, but Guseinov replaced the complete and orthonormal functions by nonorthogonal 
Slater-type functions and rearranged the resulting expansions. This is a dangerous operation whose validity must be 
checked. It is shown that the one-center limit r' = of Guseinov's rearranged Yukawa addition theorems as well as of 
several other addition theorems does not exist. Moreover, the Coulomb potential does not belong to any of the Hilbert 
spaces implicitly used by Guseinov. Accordingly, one-range addition theorems for the Coulomb potential diverge in 
the mean. Instead, these one-range addition theorems have to interpreted as expansions of generalized functions in 
the sense of Schwartz that converge weakly in suitable functionals. 

Keywords: Coulomb potential, addition theorems, exponentially decaying functions, orthogonal and nonorthogonal 
fexpansions, Hilbert spaces, generalized functions, weak convergence. 
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1 Introduction 

In many subfields of physics and physical chemistry, an 
essential step towards a solution of a problem consists in 
a separation of variables. Principal mathematical tools, 
that can accomplish such a separation of variables, are so- 
called addition theorems. These are expansions of a given 
function /(r±r') with r, r' £ M.^ in products of other func- 
tions that only depend on either r or r'. A review of the 
relevant literature with an emphasis on electronic struc- 
ture calculations can be found in [H El E] . Applications 
of addition theorems in classical physics are described in 
books by Martin [i] and Jones [S] . 

In atomic or molecular calculations, we are predomi- 
nantly interested in irreducible spherical tensors 

Frir) ^ ,U{r)Yrir/r), (1.1) 

that can be represented as products of a radial function 
fe{r) multiplied by a surface spherical harmonic Yf^{6, (p) 
defined by (|A.1|) . The convenient orthonormality and com- 
pleteness properties of spherical harmonics make it highly 
desirable that the functions of either r or r', which occur 
in the expansion of /(r ± r'), arc also irreducible spheri- 
cal tensors of the type of p.ip . Thus, addition theorems 
are expansions in terms of products of spherical harmonics 
with arguments 9,(p — r/r and 6*', ~ r' /r' , respectively. 

The best known example of such an addition theorem 
is the Laplace expansion of the Coulomb potential: 

r< I = min(r, r') , |r> | = max(r, r') . (1.2) 

Here, and Z™ denote regular and irregular solid har- 
monics defined by (|A.2|) and (|A.3p . respectively. 

The Laplace expansion ()1.2p leads to a separation of 
the variables r and r' and decouples integration variables 
in expectation values of the Coulomb potential. How- 
ever, the right-hand side of (|1.2p depends on r and r' 
only indirectly via the two vectors r< and r> that satisfy 
k<l < l''>l- Hence, the Laplace expansion has a two- 
range form, depending on the relative length of r and r' . 
This is a complication, which occurs frequently among ad- 
dition theorems and which can lead to nontrivial technical 
problems. 

As discussed in more details in [3 HI 13] , an addition 
theorem for /(r±r') can be derived by performing a three- 
dimensional Taylor expansion (see for example [SI p. 181]): 

f{r±r') = r^/W = e±'--^/(r). (1.3) 

n=0 

Thus, the translation operator 

gr'-V _ ^x'd/dx ^y'd/dy ^z'd/dz (14) 



generates f(r+r') by performing a three-dimensional Tay- 
lor expansion of / around r with shift vector r'. Since the 
variables r and r' are separated, the series expansion (|1.3p 
is indeed an addition theorem. 

We could also expand / around r' and use r as the 
shift vector. This would produce an addition theorem for 
/(r ± r') in which the roles of r and r' are interchanged. 
Both approaches are mathematically legitimate and equiv- 
alent if / is analytic in the sense of complex analysis at r, 
r', and r ± r' for essentially arbitrary vectors r,r' £ M^. 
Unfortunately, this is normally not true. Most functions 
of interest in atomic and molecular electronic structure 
theory are either singular or not analytic at the origin 
r — (0,0,0). Obvious examples are the Coulomb poten- 
tial, which is singular at the origin, or the Is hydrogen 
eigenfunction, which possesses a cusp at the origin. In 
fact, all exponentially decaying function sets discussed in 
this article are not analytic at the origin. 

Nonanalyticity is not a mathematical sophistication 
that can safely be ignored in practical applications. The 
radius of convergence of the Taylor expansion of a func- 
tion is determined by the location of its nearest singularity 
or pole. If a function / is not analytic at the origin, an 
addition theorem derived via (|1.3p can only converge if 
\r\ > |r'|, which gives rise to a two-range form. 

A two-range form can be quite inconvenient. The use 
of such an addition theorem in a multicenter integral ul- 
timately leads to indefinite integrals of special functions, 
whose efficient and reliable evaluations can be very difhcult 
(I struggled with these problems long ago in my diploma 
thesis [7 , which was published in condensed form in [Sj). 
In contrast, a Is Gaussian exp(— /3r^) with /3 > is ana- 
lytic for all r e M'^. This implies that its addition theorem, 
which can be obtained via (|1.3p , has a one-range form (see 
for example [9, Eq. (9)] or also 10, Eq. (17)]). 

The fact that addition theorems derived via (|1.3p are 
just rearranged three-dimensional Taylor series automati- 
cally implies that they converge pointwise. 

It is one of the central results of modern mathematics 
that convergence is not an abstract or global property. It 
depends very much on how we measure it. Pointwise con- 
vergence is a very demanding requirement. Moreover, it is 
not really necessary if we only want to use addition theo- 
rems in multicenter integrals. Therefore, it makes sense to 
wonder, whether computational benefits can be achieved 
by relaxing the requirement of pointwise convergence and 
by demanding instead a weaker form of convergence. 

This is indeed possible. Often, one-range addition the- 
orems can be constructed that are essentially expansions 
of /(r±r') in terms of complete orthonormal function sets. 
Such an addition theorem converges in the mean with re- 
spect to the norm of the corresponding Hilbert space. This 
had already been done by several authors (see for example 
[ni[I3[I3IIllI3[Ill[nilIl[Iljand references therein) . 
Guseinov also derived and applied one-range addition the- 
orems in several articles 20^ )21i (22j . 
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In recent years, Guseinov and his coworkers made an 
industry out of one-range addition theorems of exponen- 
tially decaying functions. They published - in addition to 
several other articles on different topics - a very long list 
of articles on the derivation and application of one-range 
addition theorems [23llM[25l[Mll23|28l|29l[30l[3ll|3a[33 

'SU. "37, 'W, W.. 'm, "ir, 'W, 1^, IB', TP, "W, 

[Ml IB71 HH [H3 [HI] , that are often remarkably similar and 
that do not give due credit to the previous work of others. 

At first sight, this long list of articles looks quite im- 
pressive. Unfortunately, from a mathematical point of 
view Guseinov's treatment of one-range addition theorems 
is at best questionable and in some cases fundamentally 
flawed. Guseinov and his coworkers ignored basic facts 
of Hilbert space and approximation theory as well as all 
questions of convergence and existence. Such a cavalier 
attitude is not uncommon among scientists, but it is com- 
pletely unacceptable in the context of addition theorems 
and multicenter integrals which are essentially a mathe- 
matical topic. 

As discussed in this article, Guseinov's approach leads 
to serious problems, which are particularly evident in the 
case of Guseinov's derivation of a one-range addition the- 
orem for the Coulomb potential [351. However, as dis- 
cussed in more details later, serious mathematical prob- 
lems persist also in all other articles dealing with one- 
range addition theorem. Therefore, the title of this ar- 
ticle is somewhat misleading. In this article, I will try to 
provide a reasonably balanced and detailed discussion of 
one-range addition theorems. Accordingly, my criticism is 
not limited to ^39J and extends to all articles by Guseinov 
and coworkers dealing with one-range addition theorems 
P?, 'M', TP, "W, '3D', '3r, T?, '33', '37, '37, '3B', "37, "SS", 
[39, 40, 41, 42, 43,, 44, 45. 46. 47. 48, 49, 50, 51, 52, 53, 54, 

[53 [SI [Ml [SS [Ml [M [Ml [M [13 [Ml [13 [Ml [13 [Ml • 

Guseinov's starting point for his one-range addition 
theorems of the Coulomb potential was a class of one-range 
addition theorems [23l [241 123 126] for Slater- type functions 
which in unnormalized form can be expressed as follows: 

X%[LiP,r) = (/3r)^-^-ic-'^'-3^f (/3r). (1.5) 

Here, N is some kind of generalized principal quantum 
number, L and M are the usual angular momentum quan- 
tum numbers, and is a positive scaling parameter. In 
the vast majority of articles dealing with Slater-type func- 
tions, N is assumed to be a positive integer > 1, i.e., 
N E N. However, several authors - among them Guseinov 
and coworkers X lU, TP, |23 [33 [S?, '33', "iT, 33", 
[S3I56, 57,M,M,M, 61. 43,Ml[Z3[ZIll22j|73,i4,I5,£6j 
- considered also nonintegral principal quantum numbers 
N £R\N. 

Guseinov's one-range addition theorems for Slater-type 
functions with nonintegral principal quantum numbers 
[211 123 [21] contain as the special case N = L — M ^ 
addition theorems for the Yukawa potential exp(— /3r) /r 



[76] , which may be viewed to be an exponentially screened 
Coulomb potential. Guseinov derived his one-range ad- 
dition theorems for the Coulomb potential by exploiting 
the obvious relationship l/r = lim^^o exp(— /?r)/r in the 
one-range addition theorems for the Yukawa potential. 

At first sight, Guseinov's derivation of a class of one- 
range addition theorems for the Coulomb potential looks 
like a remarkable achievement. So far, the Laplace expan- 
sion p.2p has routinely been used in atomic and molecular 
electronic structure calculations, but its characteristic two- 
range form can easily lead to nontrivial technical problems. 
Obviously, a one-range addition theorem for the Coulomb 
potential should simplify subsequent integrations in mul- 
ticenter integrals substantially. 

As already remarked above, one-range addition theo- 
rems are essentially expansions of a given function /(r±r') 
belonging to a suitable Hilbert space in terms of functions 
that are complete and orthonormal in this Hilbert space. 
It is an essential feature of orthogonal expansions in gen- 
eral and of one-range addition theorems in special that 
they only converge in the mean with respect to the norm 
of the underlying Hilbert space, but not necessarily point- 
wise (see for example ^773). 

The fact, that orthogonal expansions only converge in 
the mean if the function, which is to be expanded, is an 
element of the corresponding Hilbert space, has an obvi- 
ous consequence: Since the Coulomb potential does not 
belong to any of the Hilbert spaces, which Guseinov im- 
plicitly used and which all involve an integration over the 
whole M'^, Guseinov's one-range addition theorems for the 
Coulomb potential all diverge in the mean. 

But this is not the only problem. For the deriva- 
tion of his addition theorems, Guseinov had first ex- 
panded Slater-type functions xj^^(/3,r ± r') with arbi- 
trary principal quantum numbers A'^ S R in terms of 
the functions fc^'™f(/3, f) defined by (|4.16p . yielding the 
addition theorems (|6.ip . The functions j,5'™£(/3,r) with 
k — —1, 0, 1, 2, ... , which Guseinov had introduced in [Ml 
Eq. (1)], are - as discussed in Section^]- complete and or- 
thonormal in the weighted Hilbert space L^fc(IR'^) defined 
by (14.20p . As discussed in more details in Section Gu- 
seinov's functions generalize some other, well established 
function sets also based on the generalized Laguerre poly- 
nomials. 

For some reasons, which had never had been explained 
comprehensively and which I thus do not really under- 
stand, Guseinov considered it to be advantageous to re- 
place in his addition theorems (|6.ip his complete and or- 
thonormal functions fc^'™£(/3, J") by nonorthogonal Slater- 
type functions with integral principal quantum numbers 
and to rearrange the order of summations of his expan- 
sions. In this way, Guseinov obtained the expansions (|6.8p 
of Slater-type functions with in general nonintegral prin- 
cipal quantum numbers in terms of Slater-type functions 
with integral principal quantum numbers located at a dif- 
ferent center. 
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Due to the intrinsic complexity of addition theorems, 
it is by no means easy to decide whether such a rearrange- 
ment is legitimate and whether the resulting expansions 
(j6.8p are mathematically meaningful. In general, this is 
an open question, but if the principal quantum number of 
Xn LiP^ r±r') is nonintegral, N G M\N, an afhrmative an- 
swer is possible: The one-center limit r' = in (|6.8p does 
not exist, which means that Guseinov's rearranged addi- 
tion theorems (|6.8p do not exist for the whole argument 
set R3 X r3. 

This conclusion is highly consequential: The rear- 
ranged addition theorems (|6.8p with N = L = M = {) 
were the starting point for Guseinov's derivation of one- 
range addition theorems for the Coulomb potential. Thus, 
Guseinov obtained his one-range addition theorems for the 
Coulomb potential via addition theorems for the Yukawa 
potential that do not exist for the whole argument set. 

Guseinov was not the first one who had derived a di- 
vergent expansion of the Coulomb potential in terms of 
a complete and orthonormal function set. Salmon, Birss, 
and Ruedenberg [75] derived a bipolar expansion of the 
Coulomb potential in terms of the Gaussian-type eigen- 
functions of a three-dimensional isotropic harmonic oscil- 
lator defined by (|4.7|) which are complete and orthonormal 
in the Hilbert space i^(R'^) of square integrable functions 
defined by (|2.8p (see for instance jX8j Section V]). However, 
Silverstone and Kay [TS] demonstrated that this bipolar 
expansion diverges. This observation was later confirmed 
by Ruedenberg and Salmon [80] . 

Because of their divergence, one might dismiss Gu- 
seinov's one-range addition theorems of the Coulomb po- 
tential to be practically useless, just as the expansion of 
Salmon, Birss, and Ruedenberg [78] had been dismissed. 
However, this would be premature and the situation is ac- 
tually more complicated but also less hopeless than it may 
appear at first sight. As discussed in Section [3 these ad- 
dition theorems can be interpreted as weakly convergent 
expansions of generalized functions, that are mathemati- 
cally meaningful and yield convergent results when used 
in suitable functionals. 

It is quite obvious that Guseinov had failed to under- 
stand the mathematical theory behind one-range addition 
theorems. Thus is bad enough. However, the referees of 
Guseinov's numerous recent articles on one-range addition 
theorems and related topics [33[ll|23[ll[13[21[li[3ni 
[FT, Tl.'l^.'l^.'Wi^^ W, "ir, 1?, I^", '31113 

M, .47, 4.8, 49, 50j [SH [521 [53, ^ ^ M, ^ ^ ^ [6Ql 
[nillM[^[Mi[Ml[SSl[S71[Si[Ml[M] apparently also failed 
to understand this theory. Therefore, I will try to give a 
compact, but hopefully comprehensive treatment of one- 
range addition theorems and of the mathematical tools, 
which are needed for their derivation, interpretation, and 
application. 

Section [2| gives a compact review of the for our pur- 
poses most fundamental properties of Hilbert spaces in 
the context of approximation theory and of orthogonal 



expansions. Section [3] discusses the derivation and the 
basic properties of one-range addition theorems. In Sec- 
tion [H complete and orthonormal Laguerre-type function 
sets and their corresponding Hilbert spaces are discussed. 
Section [5| discusses the derivation of one-range addition 
theorems for Slater-type functions in terms of Guseinov's 
functions fe^™£(/3, r) using techniques that were applied 
by Filter and Steinborn J2J for the derivation of related 
addition theorems. Section [Hj discusses Guseinov's highly 
questionable derivation of one-range addition theorems for 
Slater-type functions. Section [7] discusses how divergent 
one-range addition theorems of the Coulomb potential can 
be interpreted as weakly convergent expansions that are 
mathematically meaningful in suitable functionals. Sec- 
tion [HI discusses alternatives to the differentiation tech- 
niques, which Guseinov had used for the generation of new 
one-range addition theorems and which are neither conve- 
nient from a technical point of view nor mathematically 
justified. This article is concluded by a short summary in 
Section [H Finally, there is Appendix \K\ listing some basic 
facts about spherical harmonics and Gaunt coefficients. 

2 Basic Hilbert Space Theory 

Let us assume that V is a vector space over the complex 
numbers C that possesses an inner product (-j-) : V x V ^ 
C (see for instance [811 p. 36]). If V is complete with re- 
spect to the norm || • || : V ^ M+ defined by 

\\u\\ = (2.1) 

i.e., if every Cauchy sequence in V converges with respect 
to (|2.ip to an element of V, then V is called a Hilbert space. 
In the mathematical literature, convergence in the means 
with respect to the norm || • || is frequently called strong 
convergence, in contrast to weak convergence which will 
be discussed later in Section [7| 

Hilbert spaces play a major role in various branches of 
mathematics and mathematical physics and in particular 
also in approximation theory. If / is an element of some 
Hilbert space Ti. and if {'/'nj^o linearly independent 
and complete in 7i, then we can construct approximations 

N 

In = (2.2) 

n=0 

to /, where is a finite integer. The coefficients C,^^ 
are chosen in such a way that the mean square deviation 
11/ - /wP = (/ - /ivi/ - In) becomes minimal. 

It is one of the central results of approximation the- 
ory that 11/ — /atIP becomes minimal if {</'n}J5°,o 
orthonormal function set satisfying {(fn\'^n') = Snn' for all 
indices n and n' and if the coefficients are chosen accord- 
ing to = {(pn\f) (see for example |S21 Theorem 9 on 
p. 51]). Since the coefficients {(pn\f) do not depend on the 
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truncation order N in fj^, f possesses an expansion 



(2.3) 



n=0 



in terms of the orthogonal function set {'Pn}?^=o that con- 
verges in the mean with respect to the norm || • || of the 
Hilbert space H. 

This result emphasizes the central role of orthogonal 
expansions both in Hilbert spaces as well as in approxima- 
tion theory. If the functions {(pn}^=o only normalized 
satisfying {ipn\^n) = 1 but not orthogonal, then it is in 
general only possible to construct finite approximations 
/at of the type of ()2.2|) by minimizing the mean square 
(/ - /wl/ - /jv), but we cannot 

(N) 



deviation \\f - /nW 

tacitly assume that the coefficients Cn" ' in Jn converge to 

oo. Accordingly 



a well defined limit C, 
the existence of an expansion 



as N 



(2.4) 



is not guaranteed if the functions {</5n}$^Lo ^ot orthog- 
onal. This fact is well documented both in the mathemati- 
cal literature (see for example [82^, Theorem 10 on p. 54] or 
[83l Section 1.4]) as well as in the literature on electronic 
structure calculations [Si l [551 [551 [57 1 [55 1 [55] ) . 

Thus, completeness of a function set {ipn}^=a in a 
Hilbert space H does not suffice to guarantee the existence 
of formal expansions in terms of these functions: These ex- 
pansions may or they may not exist. Horrifying examples 
of pathologies of nonorthogonal expansions can be found 
in [Si Section HI.I]. 

There is another, practically very consequential aspect 
of orthogonal expansions of the type of (|2.3p : The coeffi- 
cients (ipnlf) of orthogonal expansions satisfy Parseval's 
equality (see for example [5TJ Eq. (II. 2) on p. 45]) 



(2.5) 



Thus, the coefficients {ipn\f) are bounded in magnitude 
and they have to vanish as n — s- cx). This may well be the 
main reason why orthogonal expansions tend to be com- 
putationally well behaved. In the case of nonorthogonal 
expansions of the type of (|2.4|) it can instead happen that 
the expansion coefficients C„ are unbounded, have alter- 
nating signs, and increase in magnitude with increasing 
index (see for example [11, Appendix E on pp. 162 - 164] 
or [551 Table I on p. 166]). Such a behavior can easily lead 
to a cancellation of significant digits and to a catastrophic 
accumulation of rounding errors. 

If we remove a single function from a complete and 
orthonormal function set, it becomes incomplete. In the 
case of nonorthogonal expansions of the type of (j2.4p , the 
situation is much more complicated since nonorthogonal 



function sets are in general overcomplete as well as almost 
linearly dependent (see for example i87i i88) and references 
therein). Again, this can be the source of serious numerical 
problems. 

Of course, there are situations in which nonorthogonal 
expansions offer computational advantages (see fore exam- 
ple the discussion in '90]). However, in the vast majority 
of all cases, orthogonal expansions are clearly superior. 
Consequently, one should should not voluntarily abandon 
the highly useful feature of orthogonality unless there are 
truly compelling reasons. 

The natural inner product for effective one-particle 
wave functions f,g: R"^ — > C in atomic and molecular elec- 
tronic structure calculations on the basis of the Hartree- 
Fock-Roothaan equations [511 E21 US] is given by the inte- 
gral 

(/|5)2 = J [f{r)Yg{r)d'r. (2.6) 

As always in this article, integration extends over the 
whole M.^. The corresponding norm satisfies 



II/II2 = VWh- 

Accordingly, the Hilbert space 



(2.7) 



C| II/II2 <cx3} 



(2.8) 



of square-integrable one-particle wave functions provides 
the natural setting for atomic and molecular electronic 
structure calculations in general as well as for the deriva- 
tion of one-range addition theorems in special. 

The whole formalism of inner products, norms, and 
function spaces can be generalized to include weight func- 
tions w: M"^ ^ K+. If w(r) > is such a positive weight 
function, we can define an inner product with respect to 
the weight function w for functions /, g : R'^ — » C according 
to 

(/l5)u,,2 = J [f{r)Yw{r)g{r)d'r. (2.9) 



On the basis of the inner product (|2.9p . the norm of a 
function / : R'^ ^ C with respect to the weight function w 
is defined according to 



D,2 — \J (/!/)«), 2 : 



(2.10) 



and the Hilbert space L^(R'^) of square integrable func- 
tions with respect to the weight function w is defined via 
the norm (|2.10p according to 



J w{r) |/(r)p d^r < ooj 
/ll^,2<oo}. (2.11) 



By augmenting the inner product (|2.6p by a nontriv- 
ial weight function w{r) ^ 1 according to (|2.9p . it should 
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at least in principle be possible to accomplish some fine- 
tuning in approximation processes. So, it looks like an 
obvious idea to use instead of L^^jjs^ weighted Hilbert 
(M^'j ^]^g^^ jg better adapted to the problem under con- 
sideration. 

Unfortunately, there are some principal problems 
which must not be ignored. In general, we neither have 
Ll{R^) C L^(R^) nor L^{R^) C L^(M^). Thus, the 
two Hilbert spaces L'^{M.'^) and L^(R^) are mathemati- 
cally inequivalent. Accordingly, the expansion of a func- 
tion / G i2(^R3-j j-^ terms of a function set, that is com- 
plete and orthonormal in L^(R^), does not necessarily lead 
to a convergent result. Moreover, expansions of a given 
function in terms of function sets, that are complete and 
orthonormal in either L'^{M.^) or in L^(M^), respectively, 
can also be computationally different, since they can have 
substantially different rates of convergence (see the simple 
example discussed in Section [3]). 

Even the Hilbert space L'^{M.^) of square integrable 
functions is - loosely speaking - too large for electronic 
structure calculations on the basis of the Hartree-Fock- 
Roothaan equations. All effective one-particle wave func- 
tions must belong to L^^jjS-j^ j^^^ -j-j^g converse is not nec- 
essarily true, i.e., there are elements of ^^(rS) that cannot 
be used in electronic structure calculations. An example is 
Yukawa potential exp(— /3r)/r. It belongs to L'^{M.'^), but 
it cannot be used as a trial function in electronic structure 
calculations since the expectation values of the kinetic en- 
ergy and of the Coulomb potential do not exist. 

The necessary exclusion of unsuitable elements from 
£2(1^3) can be accomplished via some other generaliza- 
tions of spaces, the so-called Sobolev spaces (see for 
example |94j [95l [96]), which are of considerable impor- 
tance not only in electronic structure theory. The in the 
context of electronic structure calculations most important 



Sobolev space W2 



(1) 



can be defined via the following 



inner product for differentiable functions /, g : 



(/l5)2,i - / [/(r)] 



1* V 



C: 



(2.12) 



Here, V — [d/dx,d/dy,d/dz) is the three-dimensional 
gradient operator, and 77 is a real constant. As usual, in- 
tegration extends over the whole K.'^. 

On the basis of the inner product (|2.12p , the following 
Sobolev-type norm of a differentiable function / : M"^ ^ C 
can be defined: 



((/))2,1 = •\/(/|/)2,l- 



(2.13) 



Both the inner product {f\g)2,i as well as the norm ((/))2,i 
depend on the real parameter rj. Since it can be shown that 
the norms defined via (|2.13p are equivalent for all 77 £ M \ 
and thus give rise to the same topology, their dependence 
on r] is not explicitly indicated. 

The Sobolev space W^^\r^) is defined via the norm 



(|2.13p according to 



(1 



'(R3) 

{/:» 
{/:! 



[fir)]' 



27^2 



■fir 



d^r < 00 



C| ((/))2,1 < ^} 



} 

2.14) 



Obviously, the Sobolev space W. 



(1) 



is also a Hilbert 



.(1) 



is also a proper subspace of 



space. In addition, W2 
the Hilbert space L'^{m.^), i.e., W^^\m.^) C ^2(^3)^ Never- 



theless, approximation processes in L' 



and Wr 



(1) 



respectively, can differ substantially. Loosely speaking, we 
can say that an approximation scheme in 



imates / in the mean, but in W2^^^(]R'^) both / as well 
as V/ are approximated in the mean. Thus, in Sobolev 
spaces certain possible pathologies can be avoided. This 
can have far-reaching consequences in quantum mechan- 
ical calculations which are essentially special approxima- 
tion processes. In [85l Section 9], it was shown that com- 
pleteness of a one- particle basis in L'^{M.^) does not suffice 
to guarantee the convergence of computations based on 
the Rayleigh-Ritz variational principle. This is only guar- 
anteed if the one-particle basis is complete in W^2^''(l]^^) 
(see also 071 [SH])- 

3 One-Range Addition Theorems 



approx- 



Let us assume that / S L'^{R'^) and that the functions 
{v^eir)}n,e,m are complete and orthonormal in L'^{M.'^). 
It makes sense to assume that the functions {'P™iir)}n/,m 
are irreducible spherical tensors of the type of ()l.ip . Thus, 
the index n can be viewed to be some kind of generalized 
principal quantum number, and £ and m are the usual 
angular momentum quantum numbers. 

An addition theorem for f{r ± r'), which converges 
in the mean with respect to the norm of i2(]]j3-j^ 
constructed by expanding / in terms of the orthonormal 
functions Wn.eir)}n,£,7n- 



(3.1a) 



f{r±r') = C:],if;±r')^l],{r), 

n^rn 

CZ,if:±r') = J [vZiir)Tfir±r')d'r. (3.1b) 

The summation limit in (j3.ip depend on the exact defini- 
tion of the function set {lys™^ (?")}«/, m- This article always 
uses the convention 



(3.2) 



n=l 1=0 m=- 



which is in agreement with the usual convention for the 
bound state hydrogen eigenfunctions. 

The expansion p.ip is a one-range addition theorem 
since the variables r and r' are completely separated: The 
dependence on r is entirely contained in the functions 
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(/5™£(r), whereas r' occurs only in the expansion coeffi- 
cients Cl^g{f; ±r') which are overlap integrals. 

If the overlap integrals C''^^{f;±r') can be expanded 
in terms of the functions fl"i{r') according to 

C™,(/;±r') ^ T-,\T^,{f;±)v^,[Ar'), (3.3a) 

T:fz^'{f;±) - / (r')]*c,'::,(/;±od^r', 

(3.3b) 

then the addition theorem p.ip assumes a completely sym- 
metrical form: 

f{r±r') = ^ T„"/7„,(/;±)^™,(r)^:?;,,(r'). (3.4) 

n' i'm' 

From a purely formal point of view, the derivation 
of one-range addition theorems is a triviality. The chal- 
lenging part is the construction of computationally conve- 
nient mathematical expressions for the overlap integrals 
C^^(/;±r') in ^ or the coefficients T^/^VJ^, (/; ±) in 
In realistic applications, we cannot not tacitly as- 
sume that the use of one-range addition theorems in a 
multicenter integral necessarily leads to rapidly conver- 
gent expansions. Therefore, we must be able to compute 
C;j;^(/;±r') and/or T^/^TJ^, (/; ±) efficiently and reliably 
even for possibly very large indices (see for instance the 
convergence rates reported in (99j). 

Symmetrical one-range addition theorems of the kind 
of p.4p . which converge in the mean with respect to the 
norm of the Hilbert space L^(M'^), were constructed by 
Fiher and Steinborn [HI Eqs. (5.11) and (5.12)] and later 
applied by Kranz and Steinborn |100j and by Trivedi and 
Steinborn [99j . 

As already remarked above, a nontrivial weight func- 
tion w{r) 7^ 1 in an inner product can give more weight 
to those regions of space in which / is large, while deem- 
phasizing the contribution from those regions in which / is 
small. Accordingly, the inclusion of a suitable weight func- 
tion w can improve convergence. It is thus an in principle 
obvious idea to construct one-range addition theorems that 
converge with respect to the norm of a suitable weighted 
Hilbert space Ll{R^) defined in (pUj) . Let us therefore 
assume that / belongs to L^(R'^). Then we can construct 
a one-range addition theorem by expanding /(r ± r') with 
respect to a function set {i^Ze{r)}n,e,m that is complete 
in L^(M'^) and orthonormal with respect to the modified 
inner product (j2.9p : 

f{r±r') 

= E Tl'^m'iLy^;±)i^nAr)i;::,[,,{r'). (3.5) 

7iim 
n i'm 

For the derivation of (|3.5p . we only have to replace (|3.1I1 



by 

fir ± r') = E ^"^(^) ' (3.6a) 

n£m 

= / [V'™,(r-)]*«i(r)/(r±r')dV, (3.6b) 

and by 

C:^,(/,i«;±r') 

= E T-J^.m'{L^;±)^n',Ar'), (3.7a) 

7i'£'m' 

= y'[^:?;,,(r')]*^«(r•')C:^,(/,«;;±r')dV. (3.7b) 
As discussed in Section lU Guseinov and coworkers 

Hii sa 113 SI Hsi Ha EHi sa ia sa ED [S3 

l57l [59l [60l [6T1 [62] used for the derivation of one-range 
addition theorems and of related mathematical objects 
predominantly the functions k'^^ eiP , 'i') defined in (|4.16p 
that are complete and orthonormal in the weighted Hilbert 
space i^fc(R'^) defined in ()4.20p . As long as the functions, 
which are to be expanded, belong to the weighted Hilbert 
space L^fc(M^), this is a completely legitimate approach 
that leads to addition theorems of the type of (|3.5p which 
converge with respect to the norm (|4.19p of L'^^{M.'^). 

Unfortunately, the discussion in Section [H shows that 
the use of a weight function r*^ with fc > 1 does not neces- 
sarily lead to better results. Moreover, it is very difficult or 
even practically impossible to decide on the basis of sim- 
ple a priori considerations for which k the weight function 
r'^ produces best results. Thus, it is by no means clear 
how we can actually profit from this additional degree of 
freedom. 

There is also another annoying practical problem which 
can easily occur in multicenter integrals: Let us assume 
that f{r ± r') is expanded in terms of a function set 
{V'™£('")}n,f,m that is complete in a weighted Hilbert space 
L^(M^) and orthonormal with respect to the correspond- 
ing inner product (|2.9p . If we use this one-range addi- 
tion theorem in a multicenter integral containing /(r±r'), 
the weight function w{r), which makes (*")}«, or- 
thonormal, normally docs not occur there. Consequently, 
subsequent integrations involving f{r ± r') are not neces- 
sarily simplified by orthogonality. Thus, it may be more 
difficult to apply a one-range addition theorems of the type 
of (j3.5p based on a nontrivial weight function w{r) ^ 1 
than an addition theorem of the type of p.4p . It is not 
the complexity of an addition theorem, that really mat- 
ters, but the complexity of the series expansion, which we 
obtain by inserting an addition theorem into a multicenter 
integral. 
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If we compare one-range and two-range addition the- 
orems, it is obvious that one-range addition theorems 
greatly facihtate subsequent integrations. Thus, they seem 
to be clearly superior to two-range addition theorems of 
the type of the Laplace expansion ()1.2p . 

However, one-range addition theorems also have disad- 
vantages, and a balanced assessment of the relative merits 
of one-range and two-range addition theorems is by no 
means easy. Firstly, one-range addition theorems usually 
have a more complicated structure than two-range addi- 
tion theorems (typically, they contain one additional infi- 
nite summation). Secondly, one-range addition theorems 
normally do not converge pointwise, but only in the mean 
with respect to the norm of the underlying Hilbert space. 
This can be quite advantageous, since it makes it possible 
to expand functions with singularities and/or discontinu- 
ities. If, however, singularities and/or discontinuities are 
present, then a nonsmooth function has to be approxi- 
mated by smooth functions. Very often, this leads to slow 
convergence. 

The probably most severe disadvantage of one-range 
addition theorems compared to two-range addition theo- 
rems is that the approach sketched above cannot be ap- 
plied to all functions of interest. As remarked in Section [21 
inner products for effective one-particle wave functions in 
atomic and molecular electronic structure calculations all 
involve integrations over the whole M"^. Thus, an effective 
one-particle wave function / has to belong to the Hilbert 
space L'^{M.^) or also to a suitable weighted Hilbert space 
L^(K^). Many functions of considerable relevance such 
as the Coulomb potential or the irregular solid harmonic 
Z™(r) are only locally integrable, not square integrable 
with respect to an integration over the whole M.^ and thus 
do not belong to L^(IR'^) or to any other weighted Hilbert 
spaces discussed in this article. If we nevertheless want to 
formulate one-range addition theorems for these functions 
by expanding them in terms of a complete and orthonor- 
mal function set, we are faced with the problem that these 
one-range addition theorems diverge in the mean. 

4 Laguerre-Type Function Sets 
and Their Hilbert Spaces 

The surface spherical harmonics Y"''{9,(j)) are complete 
and orthonormal with respect to an integration over the 
surface of the unit sphere in (an explicit proof can 
for instance be found in [1011 Section III. 7. 6]). Since more 
complex Hilbert spaces can be constructed by forming ten- 
sor products of simpler Hilbert spaces (see for example 
jlOll Section 11.6.5]), we only have to find suitable radial 
functions that are complete and orthogonal with respect 
to an integration from to cxd (see also [861 Lemma 6 on p. 
31]). Thus, we more or less automatically arrive at func- 
tion sets based on the generalized Laguerre polynomials. 
The generalized Laguerre polynomials Li^\x) with 



Re(Q;) > — 1 are orthogonal polynomials associated with 
the integration interval [0, oo) and the weight function 
w(x) — a;"exp~^. They possess the following explicit ex- 
pressions [ini p. 240]: 



n + a\ X 
n — V ) v\ 



iFi(-n; a+l;x) . 



(4.1a) 



(4.1b) 



These expressions can be used to extend the definition 
of l'C\x) to complex values of a. The generalized La- 
guerre polynomials can also be defined via the following 
Rodrigues relationship ^102, p. 241]: 

L[^\x) = a;-" £! ^ re-"2:"+"] . (4.2) 
n\ dx" 

The generalized Laguerre polynomials satisfy for 
Re(a) > —1 and m,n e No the following orthogonality 
relationship [Ml P- 241]: 



r(a 



1) 



x'^e-^U-\x)L^:^\x)dx 

II'. 

(4.3) 

Henceforth, the condition Re(a) > — 1, which is necessary 
for the existence of this and related integrals, will always 
be tacitly assumed. 

The completeness of the generalized Laguerre polyno- 
mials in the weighted Hilbert space 



c 



x"'\f{x)\^dx <oo^ (4.4) 



is a classic result of mathematical analysis (see for example 
fET, p. 33], [TU3', pp. 349 - 351], or TM, pp. 235 - 238]). 

A different convention for Laguerre polynomials is fre- 
quently used in the quantum mechanical literature. For 
example, Bethe and Salpeter [1051 Eq. (3.5)] define asso- 
ciated Laguerre functions [L™(a:)]gg with n,m € Nq via 
the Rodrigues-type relationships 



[L„(a;)]gg 



dx^ 



[Ln(x)] 



BS ■ 



d" 

dx" 



Comparison of (|4.2p and (|4.5bp implies: 



(-1)™ 

(n + m)! 



BS ■ 



(4.5a) 
(4.5b) 

(4.6) 



The convention of Bethe and Salpeter [105j is also used in 
the books by Condon and Shortley [1061 Eqs. (6) and (9) 
on p. 115] and by Condon and Odaba§i [1071 Eq. (2) on p. 
189] as well as in numerous articles. 

In my opinion, the use of associated Laguerre functions 
[L™_)_„(a;)] gg defined by ([4.5p is not recommendable. It 
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follows from (|4.6p that these functions cannot express gen- 
eralized Laguerre polynomials L^"^ with nonintegral su- 
perscripts a. This is both artificial and unnecessary. For 
example, the eigenf unctions r) of the Hamiltonian 

^-2y2 _^2j,2 ^Yie three-dimensional isotropic harmonic 
oscillator contain generalized Laguerre polynomials with 
half-integral superscripts (see for example [HI Eq. (5.4)] 
and references therein): 



= fl3/2 



2{n~£~l)l 



1/2 



T{n+ 1/2) 



(4.7) 



These functions satisfy the orthonorniality condition (see 
for example \18i Eq- (5-5)]) 



[nil,{(3,r)]*n';/:,,{P,r)d^r = S,,,,,Su'S^^, (4.8) 



and are complete and orthonormal in L'^{M.^) (see also 86, 
pp. 36 - 37]). 

In this article, only the mathematical notation for gen- 
eralized Laguerre polynomials is used. Additional conven- 
tions used in physics were discussed by Kaijser and Smith 
[ini Footnote 1 on p. 48]. 

Guseinov and coworkers use the convention for La- 
guerre polynomials that is common in quantum mechanics 
(compare ^UJl Eq. (4) on p. 189] with [Ml Eq. (3)]), al- 
though they give the book by Gradshteyn and Rhyzhik 
|109j as their primary reference on Laguerre polynomials. 
This is misleading: Gradshteyn and Rhyzhik use the math- 
ematical definition of the generalized Laguerre polynomi- 
als (compare [109,, Eq. (8.970.1) on p. 1061] with g^])). 

In principle, every function set, that is complete and 
orthonormal in either L^(M'^) or also in a suitable weighted 
Hilbert space L^(IR^), can be used for the construction of 
one-range addition theorems. If we insist on an exponen- 
tially decaying function set that is complete and orthonor- 
mal in L^(M'^), then the following functions, which were 
introduced by Hylleraas |1101 Footnote * on p. 349] , ShuU 
and Lowdin [m], and Lowdin and Shull Eq. (46)], 
are the most natural choice: 



A™ (/?,r) = (2/3)3/2 



{n-£-iy. 



X e 



{n + £+iy. 



1/2 



(4.9) 



Here, /3 > is a scaling parameter. 

Lambda functions are orthonormal with respect to a 
integration over the whole R'^: 



[K2,(fi,r)\* K^:^,.{l3,r)d?r = 5 r^r.' 5 w 5 mm' ■ (4.10) 



Accordingly, they are complete and orthonormal in 



Closely related to both bound-state hydrogen eigen- 
functions and Lambda functions are the following func- 
tions which were already used in 1928 by Hylleraas [1131 
Eq. (25) on p. 478] and which are commonly called 
Coulomb Sturmians or simply Sturmians: 



*™,(/3,r) = (2/3)3/2 



{n-l-\)\ 



2n{n + t)\ 



3/2 



(4.11) 



If we replace in (|4.1ip /3 by Z/n, we obtain the bound state 
eigenfunctions of a hydrogenlike ion with nuclear charge 
Z . These eigenfunctions are orthonormal with respect to 
an integration over the whole ^ they are incomplete 
without the inclusion of the continuum eigenfunctions (see 
for example [T7] and references therein). 

Sturmians - or rather their three-dimensional Fourier 
transforms - occur also in Fock's treatment of the hydro- 
gen atom |114j . albeit in a somewhat disguised form [TSl 
Section VI]. 

It follows at once from (|4.3p that the Sturmians satisfy 
m Eq. (4.7)]: 

[^'::,{P,r)\*-^>^:,,{P,r)d?r = ^ 5r,n' Sw Smm' ■ 

' r ' n 

(4.12) 

Accordingly, the Sturmians are a complete orthogonal set 
in the weighted Hilbert space 



i?/.(»' 



J i l/(r)pd3r < ooj. (4.13) 



/r 



is not necessar- 



Unfortunately, the Hilbert space 

ily suited for quantum mechanical calculations since we 
neither have L^(R^) C Lj^^iM.^) nor i?/^(R3) C L^{R^). 
Thus, one might arrive at the premature conclusion that 
Sturmians are not particularly suited for our purposes. 

If, however, we combine the differential equation satis- 
fied by the Sturmians UHl Eq. (4.9)], 



2/3n 
r 



K'AP.r) = 0, 



(4.14) 



with the orthogonality relationship (|4.12p . we see that 
the Sturmians satisfy the Sobolev-type orthonormality [HI 
Eq. (4.10)] 



2/32 



vl':?,,(/3,r)d3r 



(4.15) 



Accordingly, Sturmians are complete and orthonormal in 
the Sobolev space 1^2'^' (R^. 

In [25, Eq. (1)], Guseinov introduced a fairly large class 
of complete and orthonormal functions which, if the math- 
ematical notation for generalized Laguerre polynomials is 
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used, can be expressed as follows: 

■(2/3)'=+3(n-£ 



-,1/2 



{n + £ + k + iy. 
X e-^'- 4'!|^+'>(2/3r) yn^Pr) . (4.16) 

The indices satisfy n e N, A; = -1, 0, 1, 2, . . . , ^ G No < 
n — 1, —i < m < i, and the scaling parameter (3 is pos- 
itive. Guseinov and coworkers used these functions quite 
extensively in the context of one-range addition theorems 

TT, 'W, TT, 'W?, '37, 'M, w 'wr, W 

El 42, 43, 44, 45, 46, 50, 47, 48, 49, 51, 53, 54, 55, 56, 57, 

If Guseinov had used the mathematical definition (|4.ip 
of the generalized Laguerre polynomials and not the un- 
necessarily restrictive convention (|4.5p , he probably would 
have noticed that the parameter A: = — 1, 0, 1, 2, . . . in his 
functions k'^^gifi, r) need not be integral and can be gen- 
erahzed to real values k £ [—1, oo). One only has to replace 
the factorial (n-|-^-|-fc-|-l)! in (j4.16p by the gamma function 
T{n + £ + k + 2). 

Guseinov's functions satisfy the orthonormality rela- 
tionship (compare also pS', Eq. (4)]) 



= 5nn' Su' Smm' ■ (4-17) 



Accordingly, Guseinov's functions are a complete and or- 
thonormal set in the weighted Hilbert space L^fc(K^) with 
fc = — 1,0,1,2,..., which is defined via the inner product 

(/l5).^2 = / [f{r)Yr'^g{r)d'r (4.18) 



and the norm 



according to 



llr^2 — \f{f\f)r'' 



(4.19) 



L2,(r3) ^ |/:R3^c| j r^|/(r)|2d3r <oo} 

= {/: M3 ^C| ||/||,.,2 < <^]- (4.20) 

For fc = — 1, the functions k^niiP^''') ^''^^ apart from a 
slightly different normalization factor identical to the Stur- 
mian functions defined by ()4.1ip . Thus, they are complete 
and orthonormal in the Hilbert space L^j^i^^) defined 
by (|4.13p or - when used in combination with the inner 
product (|2.12p with = (3 - complete and orthogonal in 

the Sobolev space ^2"^^!^^)' which is a proper subspace 
of the Hilbert space L^(R3). It is important to notice 



that Ll/^{R^) and W^^\m.^) are not identical, although 
the Sturmians and thus Guseinov's functions with k = —1 
are complete and orthogonal in both spaces. 



For A; = 0, the functions f.'if^ ^{(3 , r) are identical to 
the Lambda functions defined by (|4.9p . Thus, they are 
complete and orthonormal in L^(R3). 

For k = 1,2,3,..., the weighted Hilbert spaces 
L'^fc(M'^) are genuinely different from the Hilbert space 
L2(r3) neither have L'^{R^) C i^JM^) nor 

L^kO^^) C L^(M^). The weight function r'' in the inner 
product (|4.18p admits functions to L^fc(IR3) that are too 
singular at the origin to belong to L^(R3), and it excludes 
some functions belonging to L^(R3) that decay slowly like 
a fixed power as r ^ 00. 

If / G i^fc(M^), then it is guaranteed that the expan- 
sion 

fir) = kKMf) k^nAP,r), (4.21a) 
kK]mf) = I [,*™,(/3,r)]*r'=/(r-)d3r, (4.21b) 



in terms of Guseinov's functions converges in the mean 
with respect to the norm (|4.19p of i^j.(R3). If, however, 
/ G L^(R'^), the convergence of this expansion with respect 
to the norm ()4.19p of L^^iR^) is not guaranteed. 

For fixed k = -1,0,1,..., the functions fc^'™f(/3,r-) 
are by construction complete and orthogonal in the cor- 
responding Hilbert space ^^^(K^)^ but for k ^ the ex- 
pansion (|4.2ip does not necessarily converge with respect 
to (|4.19p for arbitrary functions / belonging to the for our 
purposes most important Hilbert space L(R3) of square 
integrable functions. Thus, the convergence of (|4.2ip is 
only guaranteed if we have / G L^(M^) n L^JR^). This is 
a complications which we have to take into account when- 
ever we do expansions in terms of Guseinov's functions 
,*™,(/3,r) withfc^O. 

It is by no means clear whether and for which k ^ 
a weighted Hilbert space L^fc(R3) provides a proper set- 
ting for bound state electronic structure calculations. It 
certainly makes sense to be cautious and to avoid poten- 
tial complications whenever possible. Therefore, I would 
refrain from using the functions ^.'^^liP, r) with k =/= un- 
less I see obvious computational benefits. Generalization 
for the sake of generalization is rarely a good idea. 

In addition, the rate of convergence of the expan- 
sion of a given function in terms of Guseinov's functions 
fc^™^(/3, r) may depend quite strongly on the choice of k. 
For fc > 1, the inner product ()4.18p gives less weight to 
the region close to the origin and greater weight to the re- 
gion away from the origin, which may or may not improve 
convergence. 

The numerical consequences of a variation of k can be 
studied by expanding the exponential 

g-x/3r (4^^-1/2 ^ g-:r/3ry-J)^^y^)^ X > , (4.22) 
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in terms of Guseinov's functions: 

oo n — 1 t 
n=l i=Q m=~e 



(4.23a) 



r)] * e-^^^ Y^iv/r) d^r . (4.23b) 



The orthonormality of the spherical harmonics imphes: 

n 1/2 



{2(3f+%n-iy. 



{n + k + 1)! 
This expression can be rewritten as follows: 



(4.24) 



(n-iy. 



1/2 



X I e 



_(2/3)'=+3(n + fc + l)!_ 



(4.25) 



With the help of the integral Eq. 7.414.8 on p. 850] 



L^^\t) At 



r(a + n+ l)(s - 1)" 



Re(Q!) > -1 



Re(s) > , 



(4.26) 



we obtain 



(l+.)t/2^fe+2 L^^+^){t)dt 

{n + k + l)\ 



x + l 



fc+3 



{n-l)\ 



X - 1 
x + l 



n-l 



(4.27) 



and 



fc^ri.o(^' P) 



X+l 

{n + k + 1)! 
2/3)fc+3(„_ 1)1 



fc+3 



- 1/2 


'x-l 


n-l 




X+l 





(4.28) 



This result shows that the expansion (|4.23p becomes triv- 
ial for X = 1 because then only the term with n = 1 is 
different from zero. However, it also shows that for x ^ 1 
the rate of convergence of the expansion (14.23^ decreases 
with increasing fc G N. 

This example does not imply that Guseinov's functions 
necessarily produce inferior convergence rates for higher 
values of k. Most likely, other examples can be constructed 
for which higher values of k improve convergence. How- 
ever, this example shows that the inclusion of the non- 
trivial weight function in an inner product can have 



significant numerical consequences. Since it is not at all 
clear whether the weighted Hilbert spaces L^fc(R'^) with 
fc ^ provide a proper setting for bound state electronic 
structure calculations, I cannot recommend to routinely 
use Guseinov's functions r) with k ^ for the 

construction of one-range addition theorems. 

This assessment does not necessarily apply to Gu- 
seinov's functions fc^'™^(/3, ''") with k = —1, which are 
up to a slightly different normalization factor identical to 
the Sturmians defined by (|4.1ip . Sturmians are complete 
and orthogonal in the weighted Hilbert space L^y^(R'^), 
whose usefulness in electronic structure calculations is 



not clear since we neither have 



nor 

L^yj,(R3) C ^^(M^). But Sturmians are also complete 

and orthonormal in the Sobolev pace W^^\r^ defined by 
(|2.14p . which is a proper subspace of L'^{M.'^). As discussed 
in Section [51 a basis set has to be complete in W2^''(I^^ 
to guarantee the convergence of variational calculations. 
Completeness in ^^(K^) does not suffice. 

Let me emphasize once more that Laguerre expansions 
converge in general only in the mean, but not necessarily 
pointwise (see for example [77] )■ Additional conditions, 
which a function has to satisfy in order to guarantee that 
its Laguerre expansion converges pointwise, were discussed 
by Szego (115, Theorem 9.1.5 on p. 246] (see also [121 Ap- 
pendix] ) . 



5 Addition Theorems for Expo- 
nentially Decaying Functions 

Pointwise convergent two-range addition theorems for all 
the commonly occurring exponentially decaying functions 
can be constructed via the two-range addition theorems 
[H I116j of the so-called B functions that were introduced 
by Fiher and Steinborn [IIll Eq. (2.14)]: 

B™,(/3,r) = [2-+Hn + iy.]-'K_,/,i(3r)yy^{Pr), 

/3>0, neZ. (5.1) 

Here, fc„_i/2 is a reduced Bessel function of half-integral 
order n - 1/2 defined by [TTHl Eqs. (3.1) and (3.2)] 



K{z) - (2/^)1/2 z-if,(z) 



(5.2) 



and Ki,{z) is a modified Bessel function of the second kind 
^102i p. 66]. If the order of a reduced Bessel function is 
half-integral, it can be expressed as an exponential multi- 
plied by a terminating confluent hypergeometric series iFi 
(see for example [HI Eq. (3.7)]). 

B functions are comparatively complicated mathemat- 
ical objects, and it is not obvious why they should offer 
any advantages compared to other exponentially decay- 
ing function sets. However, B functions possess a three- 
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dimensional Fourier transform of remarkable simplicity: 



)'''p^^^5piTTyr(-ip). (5.3) 



This is the most consequential and also the most often 
cited resuh of my PhD thesis Eq. (7.1-6) on p. 160]. 
Later, ()5.3p was pubhshed in [1211 Eq. (3.7)]. Indepen- 
dently and almost simultaneously, (|5.3p was also derived 
by Niukkanen [T^, Eqs. (57) - (58)]. 

The exceptionally simple Fourier transform (|5.3p gives 
B functions a unique position among exponentially decay- 
ing functions. It also explains why other exponentially 
decaying functions like Slater-type functions with integral 
principal quantum numbers, bound state hydrogen eigen- 
functions, and other functions based on generalized La- 
guerre polynomials such as Lambda functions, Sturmians, 
or Guseinov's functions can all be expressed in terms of fi- 
nite linear combinations of B functions (details and further 
references can be found in [18, Section IV] or [2^ Section 
4]). Two-range addition theorems for the exponentially 
decaying functions mentioned above can be written down 
immediately by forming finite linear combinations of the 
corresponding addition theorems of B functions [21 1116j . 

B functions also turned out to be extremely useful for 
the derivation of one-range addition theorems. Filter and 
Steinborn [TU Eqs. (5.11) and (5.12)] derived symmetri- 
cal one-range addition theorems of the type of p.4p for 
Lambda and B functions by expanding them in terms of 
Lambda functions: 



nin2 

ili2 mi 



(5.4) 



(5.5) 



These addition theorems contain Gaunt coefficients de- 
fined by (|X6|l . Filter and Steinborn [12, Eqs. (4.9), (4.10), 
and (4.33)] were able to derive explicit expressions for the 
coefiicients T^,^g'^^^ and a^llf^^^ occurring in ^ and 
(|5.5p . respectively. In [TSl Eq. (7.8)] it was shown that the 
explicit expression for the coefficients T^^^^'^'^ in the addi- 
tion theorem (|5.4p can also be derived via the weakly con- 
vergent expansion of the plane wave in terms of Lambda 
functions [HI Eq. (4.38)1. The same approach works also 
for the coefficients a^^^^ """^^ in the B function addition 
theorem (|5.5p (see UM Eq. (7.7)]), and it should also work 
in the case of Slater-type functions with integral princi- 
pal quantum numbers defined by p.5p . although this has 
apparently not been done yet. 



The derivation of the addition theorems (|5.4p and (|5.5p 
by Filter and Steinborn was based on their remarkably 
compact convolution theorem of B functions [1231 Eq. 
(4.1)] which, however, can be derived more easily with the 
help of the Fourier transform ([5.3p [OH Section V]: 



/ BZ\iM [r^r'])BZ:,,iP,r')d'r' 



47r 



•^^^ {£mi + m2\£imi\i2m2) 



Ae 

E 

t=0 



(-1)* 



Ai 



V r>mi+m2 



Ui+i2-e-t+iA^^''') ■ (5-6) 

The abbreviation A£ is defined by (|A.9ll . For an appli- 
cation of this highly convenient expression for the con- 
struction of the addition theorem l|5.4p . Filter and Stein- 
born only had to express Lambda function in terms of B 
functions. This can be accomplished with the help of the 
following expression [120|, Eq. (3.3-35)] (see also [H Eq. 
(3.17) and Ref. [23] on p. 2736]): 



4")(2z) = {2n + a + l) 



i^2yTin + a + u + 1) 
^ I'Un - iy)ir(a + 2iy + 2) 



E 



Wi/2(^)- (5.7) 



This yields a finite sum representation of Lambda func- 
tions in terms of B functions 12. Eq. (3.18)]: 



A™,(/3,r-) = (2/3)3/2 2^ gilii 



{n + £+iy. 



D! 



1/2 



(^" + ^+l)>^ in + e + 2)^ 

!/=0 



u\ {£ + 5/2), 



{n^£- 



(5.8) 



In this way, all overlap integrals of Lambda functions or 
between B and Lambda functions, which occur according 
to p.ip as intermediate steps in the derivation of the one- 
range addition theorems (|5.4p and ([5.5p . can because of 
(|5.6p be expressed in terms of B functions. 

Now, one only has to express all B functions in terms 
of Lambda functions. This can be done with the help of 
the following relationship [124', Eq. (B.2) on p. 214]: 



fc„+i/2(^) - (-2)-"n!e--4-2"-i)(2z). 



(5.9) 



As is well known, a generalized Laguerre polynomial with 
superscript /3 can be expressed as a finite sum of gener- 
alized Laguerre polynomials with superscript a \W2\ p. 
249]: 



Lif\-) - E 



(/? 



-(") 



(5.10) 



This yields the following finite sum representation of a re- 
duced Bessel function with half integral order in terms of 
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generalized Laguerre polynomials with an essentially arbi- 
trary superscripts a: 

nl 

kn+l/2iz) 



2" 



m— 



2n + a + l 



Lt\2z). (5.11) 



Filter and Steinborn [T?, Eq. (3.13)] gave an expression 
which wrongly contains the additional factor (—1)". 

Thus, we can express a B function as a finite sum of 
Lambda functions (see also 12, Eq. (3.14)]): 



BZAfi,r) = (2/3)- 



^ in + 21- 



_3/2 (n + 2£ + 3)„-i 



{v + 2l + 2)\ 



3), 



-,1/2 



X A 



iP,r), nSN, /3>0. (5.12) 



Starting from the Lambda function addition theorem 
a symmetrical one-range addition theorem of the 
type of (I3.4|l for Slater-type functions with integral prin- 
cipal quantum numbers n G N can be derived easily. As 
is well known, an integral power with m € Nq can be 
expressed as a finite sum of generalized Laguerre polyno- 
mials Eq. (2) on p. 207]: 



(a + 1),; 



n=0 



(a + 1)„ 



(5.13) 



Thus, a Slater-type function with integral principal quan- 
tum number n £ N defined by (|1.5p can be expressed as a 
finite sum of Lambda functions [12 Eq. (3.10)]: 



xZiiP,r) = (2/3) 



-3/2 (2^ + 3)n-£- 



(2^ + 3), 



(iy + 2^ + 2)! 



-,1/2 



X A' 



(5.14) 



Accordingly, a symmetrical one-range addition theorem 
for Slater-type functions with integral principal quantum 
numbers n G N in terms of Lambda functions can be writ- 
ten down immediately. The expansion coefficients of this 
addition theorem are just simple finite sums of the coeffi- 



cients Tl^llf^ in (lOD . 



Alternatively, one could just as well start from the B 
function addition theorem (|5.5p since a Slater-type func- 
tion with an integral principal quantum number can also 
be expressed as a finite sum of B functions |1171 Eqs. (3.3) 
and (3.4)]. 

= 2" y (-1)- 



cr>0 

X {n-a)\BZU-aAP^r). 



(5.15) 



If the principal quantum number n is a positive inte- 
ger, this expansion terminates because of the Pochhammer 
symbols (— [n — ^ — 1]/2)(t and (— [n — P\/2)cr after a finite 
number of steps. 

Thus, the expansion coefficients of the Slater addition 
theorem can also be expressed as simple finite sums of the 
coefficients all\f+^^ in (l531) . 

The finite sum formula (|5.13p can be generalized to 
nonintegral powers with /i G C \ Nq. With the help of 
[W^, Eq. (7.414.7) on p. 850] 



f 

Jo 



n!r(a-f-l) * 
X 2Fi{~n,P+\-a + \\l/s) , 

Re(/3)>-l, Re(s)>0, 



(5.16) 



we obtain: 



r(^ + a + 1 



1) ^„ 



(-/^)n 
{a+l)n 



^gC\No, Re{n + a),Re{a) > -1. (5.17) 

If we set fj, = m with m G Nq, the infinite series on 
the right-hand side terminates because of the Pochham- 
mer symbol (— m)„ and we obtain the finite sum (|5.13p . 
However, there is a fundamental difference between (|5.13p 
and (j5.17p . The finite sum formula (j5.13p is a relationship 
among polynomials. Therefore, it is valid pointwise for ar- 
bitrary a; G C. In the case of the infinite series ()5.17p . we 
know that it converges in the mean in the radial Hilbert 
space (|4.4|) , but we cannot assume that it converges point- 
wise. Moreover, the index dependence of the series coef- 
ficients on the right-hand side of (|5.17p indicates that the 
convergence of this expansion is slow if it does not termi- 
nate. These convergence problems can be demonstrated 
by considering the special case fi = —I and x — » in 
(|5.17p . Then, the left-hand side approaches +oo. For an 
analysis of the behavior of the right-hand side of (|5.17p as 
a; ^ 0, we use (|4.1b|) and obtain 



(5.18) 



Inserting this into (j4.1bp yields: 



lim — 

x^O X 



(1)„ {a + l)n 



a 



iy ilk. 

n=0 



(5.19) 



Since (1)„ = n!, the series diverges for a; ^ to -t-cxo, 
although each individual term of the series (|5.17p with 
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/z = — 1 is well behaved as a; — > 0. Of course, pointwise 
convergence - or also divergence - of an expansion like 
(|5.17p for nonintegral /i G M \ No is not really important if 
we want to use it in integrals. However, the at best slow de- 
cay of the terms on the right-hand side of ()5.17p indicates 
that expansions of a Slater-type function with nonintegral 
principal quantum number n ^ N in terms of Lambda 
functions or other exponentially decaying Laguerre-type 
function sets converge slowly. 



rived via (|5.10p : 



By combining (|5.17p with the Lambda function addi- 
tion theorem (|5.4p , a symmetrical one-range addition the- 
orem for a Slater-type functions with a nonintegral prin- 
cipal quantum numbers in terms of Lambda functions can 
be formulated easily. The only principal constraint on the 
nonintegral principal quantum number is that the corre- 
sponding Slater-type function must be square integrable, 
i.e., it must belong to the Hilbert space L^(R^). How- 
ever, the expansion coefficients of this one range addition 
theorem are now given by an infinite series involving the 
coefficients T^^g^^^ in ((5Ti)) . As remarked above, we have 
no a priori reason to assume that the convergence of this 
series would be rapid. So, unless somebody accomplishes 
a substantial simplification - for example by expressing 
an inner sum in closed form via a summation theorem of 
a generalized hypergeometric series with unit argument - 
we are confronted with a one-range addition theorem for 
Slater-type functions with nonintegral principal quantum 
numbers that is significantly more complicated and less 
suited for practical work than the corresponding addition 
theorem for Slater-type functions with integral principal 
quantum numbers. 

Let us now assume that we have a symmetrical one- 
range addition theorem for Slater-type functions in terms 
of Lambda functions of the type of (15. 4p or (|5.5p . Then 
we only have to replace the Lambda functions by Gu- 
seinov's functions to obtain an expansion of a Slater-type 
function in terms of Guseinov's functions r) with 

k = -1,0,1,2,... defined by ^J^. With'the help of 
(|5.10p we obtain the following expression for a Lambda 
function as a finite sum of Guseinov's functions: 



min(n-£-l,fc) x -il/2 



{n + l + 2)k- 
r). (5.20) 



In this way, expansions in terms of Lambda functions can 
be transformed easily to expansions in terms of Guseinov's 
functions. 

An inverse relationship - the expansion of Guseinov's 
functions in terms of Lambda functions - can also be de- 



n-t-l 



,*™,(/3,r) = (2/3)^-/2 ^ 



1/2 



(5.21) 



With the help of this finite sum, we could - starting from 
the addition theorem (|5.4p for Lambda functions - im- 
mediately write down a symmetrical one-range addition 
theorem for Guseinov's function. Moreover, all Lambda 
functions in this addition theorem could be replaced by 
Guseinov's functions via (|5.20p . 

The approach of Filter and Steinborn [12], which is 
based on the convolution theorem (|5.6p of B functions, 
can also be used to construct from the scratch symmet- 
rical one-range addition theorems that are expansions in 
terms of Guseinov's functions. Setting w{r) ~ r^ and 
^/;™^(r) = fc^'™£(/?,r-) in ([33]) yields the following one- 
range addition theorem: 



/(r±r') 



E 



fcT"^^m'(/; /3, ±) 



X fc*™ (/3,r),vl;™,,,(/3,r'), (5.22a) 
fcT"^™„i'(/;/3, ±) 

[,*">(/?, r')] * (r')'= ,C;:;,(/; /3, ±r') dV , 

(5.22b) 

,C™,(/;/3,±r') 

= / [XAP,r)yr'^f{r±r')d^r. (5.22c) 

If / £ L^fc(K'^), this addition theorem converges in the 
mean according to the norm (|4.18p of the weighted Hilbert 
space L2,(R3). 

For the derivation of an addition theorem of the type 
of (|5.22p for Guseinov's functions or B functions via the 
convolution theorem (|5.6p of B functions, we have to ex- 
press Guseinov's functions in terms of B functions. This 
can be done with the help of (|5.7p . yielding 



{n + e + k + 1)1 



2^+1 {n-£-iy. 
{2n + k + l)r{l/2) {£+l)\ 
r(^ + 2 + A:/2)r(^ + [fc-f 5]/2) 



1/2 



E 



{-n + £+l)^{n + e + k + 2)^ {I -)- 2)^ 



^!(£ + 2 + fc/2)^(£+[fc + 5]/2)^ 
X S::;i^,(/3,r). (5.23) 

With the help of this relationship, two-range addition the- 
orems of Guseinov's function can be written down imme- 
diately by forming finite linear combinations of the corre- 
sponding addition theorems of B functions [21 1116j . 
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The weight function r'^ in (|5.22b[l and ()5.22c|) , which is 
responsible for the orthonormahty of Guseinov's functions 
according to (|4.17p . can be absorbed with the help of [123, 
Eq. (6.1))] 



(-s/2),(-n- [s-l]/2). 



a 



a>0 

X kn+s-rT-l/2{z) , 

s = -1,0, 1,2, ... , neN, (5.24) 



yielding 



= {2/f3r J2 ("1 



{-s/2),{-n-[s-l]/2)„ 
cr! {n + e+l),_^ 



a>0 

s = -l,0,l,2,... , neN. (5.25) 



If s is an even integer, the Pochhammer symbol (— s/2)cr 
causes a truncation of the summation after a finite number 
of steps, and if s is odd, this truncation is accomplished 
by the Pochhammer symbol (— n — [s — l]/2)o-. 

Finally, we have to express B functions as finite linear 
combinations of Guseinov's functions. With the help of 
(|5.1ip . we obtain: 



(n + 2£ + fc + 3)„-i 

22n+2<!-l + 



E 



(n + 2£ + fc + 3)^ 



{v + 2l + k + 2)\ 



{2(3Y+^ v\ 



X fe*"V£+l/(A^). 



1/2 



(5.26) 



With the help of (jOg]) . (jOS]) . and (jOS]) it is pos- 
sible to generalize the approach of Filter and Stcinborn 
[T2], which was based on the convolution theorem (|5.6p 
of B functions and which produced expansions in terms 
of Lambda functions, to one-range addition theorems for 
exponentially decaying functions that are expansions in 
terms of Guseinov's functions j^^\l'™^(/3, r). 

We cannot expect that these addition theorems are 
as compact as the corresponding expansions in terms of 
Lambda functions, which can be obtained by setting fc = 
in Guseinov's functions. For example, (|5.23p simplifies 
considerably for fc = — 1 and fc = 18, Eqs. (4.19) and 
(4.20)]. Moreover, for fc = -1 and fc = 1, 2, . . . , the weight 
function has to be absorbed via (j5.24p . which produces 
an additional inner sum. So, unless some additional sim- 
plifications can be found, one-range addition theorems in 
terms of Guseinov's functions will have a more compli- 
cated structure than the corresponding addition theorems 
in terms of Lambda functions. 

This higher complexity is one reason why I am not in- 
terested in explicitly constructing symmetrical one-range 



addition theorems that are expansions in terms of Gu- 
seinov's functions with fc > 1. I doubt that these addition 
theorems would be useful enough to justify the effort. A 
second reason is that I have - as outlined in Section [1]- se- 
vere doubts whether the weighted Hilbcrt spaces ^^^(M^) 
with fc 7^ are really suited for bound state electronic 
structure calculations. 



6 Guseinov's One-Range Addi- 
tion Theorems for Slater- Type 
Functions 

As discussed in Section [Sj symmetrical one-range addition 
theorems of the type of (|3.ip for exponentially decaying 
functions can be derived comparatively easily via the re- 
markably compact convolution theorem (|5.6p of B func- 
tions. However, Guseinov preferred to proceed differently. 

In [23] , Guseinov derived one-range addition theorems 
for Slater-type functions with integral principal quantum 
numbers by expanding them in terms of Sturmians and 
Lambda functions, and in [21j he did this for Slater-type 
functions with nonintegral principal quantum numbers. In 
[25 , Guseinov introduced his functions fc^'^^(/3, r) defined 
by (|4.16p and used them for the construction of one-range 
addition theorems for Slater-type functions with integral 
principal quantum numbers. Since Sturmians and Lambda 
functions are special cases of Guseinov's functions with 
fc = — 1 and fc = 0, respectively, these addition theorems 
generalize earlier expansions in terms of Lambda functions 
and Sturmians derived in 23J. In 26] Guseinov derived 
one-range addition theorems in terms of his functions for 
so-called "central and noncentral potentials" which are 
nothing but special Slater-type functions with integral and 
nonintegral principal quantum numbers. In [27j . Guseinov 
provided a "unified treatment" of multicenter integrals of 
Slater-type functions with integral and nonintegral princi- 
pal quantum numbers by expanding Slater-type functions 
in terms of his functions. In [28|, Guseinov constructed 
addition theorems for his functions both in the coordinate 
as well as in the momentum representation by expanding 
his functions in terms of Slater-type functions with inte- 
gral principal quantum numbers. In |29j . Guseinov used 
his addition theorems for Slater-type functions to handle 
multicenter integrals of what he calls "central and noncen- 
tral interaction potentials" which are special Slater-type 
functions with integral and nonintegral principal quantum 
numbers. In [30], Guseinov provided again a "unified ana- 
lytical treatment" of multicenter integrals of "central and 
noncentral interaction potentials" via one-range addition 
theorems of Slater-type functions derived with the help of 
his functions. In [31 , Guseinov provided in this way a 
"unified analytical treatment" of multicenter nuclear at- 
traction, electric field and electric field gradient integrals 
over Slater-type functions. In [32l I35j . Guseinov provided 
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another "unified treatment" of essentially the same in- 
tegrals as in j31j . but this type he emphasized the use 
of his functions. In [33j . we find another "unified treat- 
ment" of multicenter integrals of Slater- type functions, 
but this time the potentials are called "integer and non- 
integer u Yukawa-type screened Coulomb type potentials" 
(again, these potentials are nothing but special Slater-type 
functions with integral and nonintegral principal quantum 
numbers). In [34 , Guseinov constructed one-range addi- 
tion theorems for derivatives of Slater- type functions, and 
in [361 I37j . he considered again multicenter integrals of 
"central and noncentral interaction potentials". In [38] , 
Guseinov provided another "unified treatment" , but this 
time of electronic attraction, electric field, and electric- 
field gradient multicenter integrals of "screened and non- 
screened Coulomb potentials". In [311, Guseinov derived 
one-range addition theorems for the Coulomb potential 
starting from his addition theorems for the Yukawa po- 
tential which is essentially a special Slater-type function 
Xn,l with N = L = M = 0. In Guseinov used 

the momentum space addition theorems derived in [28j for 
the construction of momentum space addition theorems 
for Slater- type functions. In [HI [l2l 031 [45l [46] , Guseinov 
considered one-range addition theorems for "Yukawa-like 
central and noncentral interaction potentials and their 
derivatives" , for derivatives of his functions, for "deriva- 
tives of integer and noninteger u Coulomb- Yukawa type 
central and noncentral potentials" , and for "combined 
Coulomb and Yukawa like central and noncentral inter- 
action potentials and their derivatives" , respectively. In 
[47] . Guseinov considered one-range expansions for two- 
center charge densities of Slater-type functions with inte- 
gral and nonintegral principal quantum numbers. Finally, 
in [35] Guseinov provided another "unified treatment" of 
expansion theorems and one-range addition theorems of 
complete orthonormal sets of functions in coordinate, mo- 
mentum and four-dimensional spaces. 

These examples show that Guseinov used quite a few 
different names for the functions he expanded, but with 
few exceptions he concentrated on one-range addition the- 
orems for Slater-type functions with integral and noninte- 
gral principal quantum numbers in the coordinate repre- 
sentation. Consequently. I will also focus on these addi- 
tion theorems and will examine critically their derivation 
by Guseinov. 

Guseinov exclusively derived unsymmetrical addition 
theorems of the type of (|3.ip or (|3.6p . in which the two 
vectors r and r' are treated differently. He never derived 
completely symmetrical addition theorems of the type of 
((331), (EUl), or ^M)- 

Moreover, Guseinov expanded Slater- type functions in 
terms of functions with a different scaling parameter. Nor- 
mally, this is a bad idea since overlap integrals with dif- 
ferent scaling parameters are significantly more compli- 
cated than overlap integrals with equal scaling parame- 
ters (this is the reason why exclusively addition theorems 



with equal scaling parameters were discussed in Section 
[5]). However, Guseinov needed this additional degree of 
freedom since he wanted to construct one-range addition 
theorems for the Coulomb potential by representing it as 
the limiting case of the Yukawa potential according to 
1/r = lim^^o exp(— /3r)/r. 

Guseinov's work on one-range addition theorems is 
based on expansions of Slater-type functions with integral 
or nonintegral principal quantum numbers in terms of his 
functions fc*"f(7'''') '^i*^ ^ ^ -1,0,1,2,... and 7 > 
(see for example [25l Eqs. (11) and (12)]): 

= E kKt^'^^'f^'^^') (6.1a) 

.X^^™"(7,/?,±r-') 

- / U*™,(7,r-)]*r'xj^^.i(/3,r-±r')dV. (6.1b) 

As remarked above, Guseinov did not derive completely 
symmetrical addition theorems of the type of ()5.22|) by 
expanding the overlap integrals (j6.1bp in terms of his func- 
tions. This is not completely satisfactory: In the Slater- 
type function xjlf ^(/3, r±r'), the vectors r and r' play an 
identical role, but not in the addition theorems (|6.ip . 

Nevertheless, a definite assessment of the relative mer- 
its of symmetrical and unsymmetrical one-range addition 
theorems is not so easy. Since we cannot tacitly as- 
sume that one-range addition theorems necessarily lead 
to rapidly convergent expansions for multicenter integrals, 
the crucial question is whether it is easier to compute ei- 
ther overlap integrals like (j6.1bp or purely numerical coef- 
ficients as they occur in (|3.3b|) . (j3.7bp . and (|5.22bp both 
effectively and reliably for large indices. 

It also makes a big difference if one has to integrate 
over only one of the two vectors r and r' or over both of 
them. For example, in [321 Table I] Guseinov considered 
three-center one-electron nuclear attraction integrals 

I{f,g;A,B,C) 

= I [/(r.-A)]*^^.9(r-C)d3r, (6.2) 

where / and g are Slater- type functions, and the vectors 

A, B, and C are atomic centers. The integral (|6.2|) is 
the most complicated one-electron integral occurring in 
molecular electronic structure calculations based on the 
Hartree-Fock-Roothaan equations. By means of a shift 
of origin, one of the three vectors A, B, and C can be 
made to vanish. Thus, two addition theorems are needed 
to decouple the arguments in (|6.2|) . 

In the case of a one-electron integral like (|6.2p . the use 
of unsymmetrical one-range addition theorems of the type 
of (|6.ip may well be feasible: It should not matter too 
much whether the remaining two of the three vectors A, 

B, and C occur in overlap integrals or in complete and 
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orthonormal functions. Nevertheless, a comparison of Gu- 
seinov's results [39, Table I] with those reported by Boufer- 
guene and Jones [126] . who had used two-range addition 
theorems, would have been of considerable interest. 

The most difficult multicenter integrals occurring in 
the LCAO-MO approach arc the notorious six-dimensional 
two-electron integrals 



^N.L.M 



(7,/3,±r') 



nr)9{r') 



d> d>' , 



(6.3) 



which describe the Coulomb interaction of two in general 
nonclassical charge distributions [/(t)]* and g{r') consist- 
ing of products of effective one-particle wave functions lo- 
cated at different atomic centers. We can only benefit in 
both integrations from the simplifying power of orthogo- 
nality if we use in (|6.3p symmetrical one-range addition 
theorems of the type of ([33), JM]), or ([5?22ll . 

In principle, this applies also to multicenter inte- 
grals of those functions which Guseinov and coworkers 
had called "central and noncentral interaction potentials" 

[^aisaEniiMiEgEaiiiiiiaiiaiisiiigE^ 

and which are nothing but special Slater-type functions 
with integral and nonintegral principal quantum numbers. 
In the case of one-electron integrals, unsymmetrical addi- 
tion theorems of the type of (|6.ip may well be sufficient, 
but if two-electron integrals of these interaction potentials 
should indeed be physically meaningful and have to be 
computed (see also |127| ). I am skeptical about the feasi- 
bility of Guseinov's approach based on the exclusive use 
of unsymmetrical addition theorems. 

The central computational problem of Guseinov's un- 
symmetrical one-range addition theorems (j6.1|) is the effi- 
cient and reliable evaluation of the overlap integrals (j6.1bp 
between Slater-type functions and the expansion functions 
j.^'™^(7,r) even for large indices. 

Guseinov's solution is simple but not necessarily good. 
It follows at once from the explicit expression (|4.ip of the 
generalized Laguerre polynomials that Guseinov's func- 
tions can be expressed as finite sums of Slater-type func- 
tions with integral principal quantum numbers: 

fe*".(7,r) = E fcG^^'(7)xi:Wi,.(7,r), (6.4a) 



{n-t-l)\ 
(-n + <'+ 1)^2"^ 



1/2 



Accordingly, the overlap integrals (|6.1bp can be expressed 
as finite sums of overlap integrals of Slater-type functions: 

XN,L,M/ a I i\ -k 
. uAra (7,/3,±r) = 7 

n-l-1 

X E feG("''^(7)S;Z;^+t+i,.,™(7,/3,±r'), (6.5a) 

v=0 



r ± r') d^' 



(6.5b) 



Using this in the addition theorems (|6.ip yields: 



X^JAP^r±r') ^ 7"" E k'fnAl^r) 

n-l-1 

X E ;.■G^'H7)S;r4;\^+l,,,„(7,A±r'). (6.6) 

i/=0 

Superficially, this approach, which is mathematically com- 
pletely legitimate, seems to have the advantage that ex- 
isting programs for overlap integrals of Slater-type func- 
tions can be used for the evaluation of the overlap integrals 
((6lb)) (this may be the reason for Guseinov's approach). 
However, stability problems are likely. The explicit expres- 
sion ()4.ip of a generalized Laguerre polynomial lI"'' (x) in 
powers of x becomes numerical unstable for larger values 
of the index n. Because of these stability problems, or- 
thogonal polynomials are normally computed recursively. 

So, as long as the opposite is not explicitly proved, it 
certainly makes sense to assume that an expression like 
(|6.4p inherits the stability problems of (|4.ip . In addition, 
conventional programs for overlap integrals of Slater-type 
functions, as they are for instance used in semiempirical 
calculations, normally cannot be used in the case of (very) 
large principal and angular momentum quantum numbers. 
Thus, we need alternative expressions for the overlap inte- 
grals (j6.1bp not based on (|6.4p . which permit an efficient 
and reliable evaluation even for (very) large indices. 

In his desire to reduce his whole formalism of one-range 
addition theorems to Slater-type functions, Guseinov even 
expressed the functions k^'^i{"i,r) on the right-hand side 
of (|6.6p by Slater-type functions according to ()6.4p (see for 
example ^ Eq. (14)]): 



7-'^E E kG^:^''\^)xT'+,+^Al,r) 

nim u' —0 
n~£~l 

X E .G(-^n7)S:);^^\^+,,,„(7,A±r'). 



(6.7) 



i/=0 



This is still correct since the sum over j^' is nothing but 
the function fc^'™£(7,»") in disguise. But Guseinov rear- 
ranged the order of summations in (|6.7p , obtaining expan- 
sions for Slater-type functions Xn hiP^''' ^ '^ith inte- 
gral or nonintegral principal quantum numbers in terms 
of Slater-type functions X^+i lif^^''') with integral princi- 
pal quantum numbers located at a different center (see for 
example Eq. (15)]): 



x'lL{P.r±r') 



E 

nim 
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We then obtain for the hypergeometric series in (|6.9p : 



n+n' -l-l 

X S:r4;t+i/.™(7,/3,±r-')- (6.8) 



This step is potentiaUy disastrous. A rearrangement of 
the order of summations of a double series is not always 
legitimate and can easily lead to a divergent result (see 
for example |128l Chapter V]). So, when going from (|6.7p 
to (|6.8p . one cannot tacitly assume that the rearranged 
expansion (|6.8p indeed converges. 

As discussed in Section a formal expansion of the 
type of ()2.3p of a function / in terms of a given func- 
tion set {(pn}'^=o exists and converges in the mean if / 
belongs to the underlying Hilbert space and if the expan- 
sion functions {<i5n}.J5Lo ^'"^ complete and orthonormal in 
this Hilbert space. If the expansion functions are only 
complete, but not orthogonal, it is possible to construct fi- 
nite approximations of the type of (|2.2p by minimizing the 
mean square deviation ||/ — JnW^ = if ^ ^ In), but 
the existence of formal expansions of the type of (|2.4p in 
terms of nonorthogonal functions is not guaranteed: These 
expansions may or may not exist. 

Slater-type functions are complete in all the Hilbert 
space considered in this article (see for example [SB, Sec- 
tion 4]), but not orthogonal. Thus, it is not clear whether 
Guseinov's rearranged addition theorems (|6.8p are mathe- 
matically meaningful. This has to be checked. 

Convergence and existence problems due to a rear- 
rangement can be illustrated via (|5.17p . which expresses 
a nonintegral power x'^ with /i e C \ Nq as an infinite 
series of generalized Laguerre polynomials and which con- 
verges in the mean in the Hilbert space (|4.4p . If we insert 
the explicit expression (|4.ip of the generalized Laguerre 
polynomials into (|5.17p and interchange the order of sum- 
mations, we obtain after some algebra: 

^ _ T{^i + a + l) 
Tia + 1) 

OO f \ k 

It looks as if we succeeded in constructing a power series 
expansion for the nonintegral power x'^ . However, the gen- 
eralized hypergeometric series i_Fo with unit argument is 
the limiting case ?/ — > — 1 of the following binomial series 
satisfying [101 p. 38] 



m=0 ^ 

(6.10) 



iFo{k-fi;l) 

{oo , /i < , 

0, k<fi>0, (6.11) 

CXD , fc > /X > . 



Thus, the power series ()6.9p does not exist and the La- 
guerre series (|5.17p for x^^ with /i e C \ No cannot be refor- 
mulated as a power series in x. Of course, this is not really 
surprising: The general power function with z G C and 
G C \ No is not analytic at z = in the sense of complex 
analysis. For /i = m with m G No, Taylor expansion of z™ 
around z = yields the trivial identity z™ = z™. 

One-range additions theorems for exponentially decay- 
ing functions are fairly complicated mathematical objects. 
Consequently, explicit proofs of their convergence or diver- 
gence are very difficult. It is understandable that Guseinov 
was not interested in presenting such proofs, although he 
should have done so in order to justify his manipulations. 
Fortunately, valuable insight can in some cases be gained 
by considering not the addition theorems themselves, but 
their much simpler one-center limits. 

Let us therefore assume that we succeeded in construct- 
ing either symmetrical addition theorems of the type of 
(j5.22p or also unsymmetrical addition theorems involving 
overlap integrals for some function /(r ± r') by expand- 
ing it in terms of Guseinov's functions. We now consider 
the one-center limit by setting r' — 0. Then, our addition 
theorem must simplify to yield an identity for /(r). Under 
fortunate circumstances, the mathematical nature of this 
identity can provide valuable insight. 

First, we set /3 = 7 and r' = in the addition theorems 
(|6.ip for Xj^ l{(3, r ± r'). After the cancellation of all com- 
mon factors we obtain expressions that are special cases of 
the expansion (|5.17p expressing an in general nonintegral 
power x'* as an infinite series of generalized Laguerre poly- 
nomials. As shown in (|6.9p - (|6.1ip . this series can only be 
reformulated as a power series in a; if is a nonnegative 
integer, fj, — m with to G No, yielding the trivial identity 
a;™ — X™. If we have instead /i G M \ No, a rearranged 
power series in x does not exist. 

Thus, we can conclude that for /3 = 7, the one-center 
limit r' = of the rearranged addition theorem (|6.8p does 
not exist if the principal quantum number TV of the Slater- 
type function Xn hif^' ^ '"'') nonintegral, A G M \ N. 

Next, we assume /3 7^ 7 in the addition theorem (|6.ip 
and set r' = 0. We then obtain after the cancellation of all 
common factors expressions that are special cases of the 
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following expansion: 



n=0 ^ ^ 

/iGK, Re(^ + a)>-l, ue(-oo,l/2). (6.12) 

This expansion, which can be derived with the help of 
(j5.16p . converges in the mean with respect to the norm of 
the weighted Hilbert space Ll_^^^{R+). For u = 0, 
simplifies to give (|5.17p . This can be shown easily with 
the help of Gauss' summation theorem !102|, p. 40]. 

If we insert the explicit expression (|4.ip of the gener- 
alized Laguerre polynomials into (|6.12p and interchange 
the order of summations, we also obtain a formal power 
series in x. Unfortunately, an analysis of the resulting 
power series becomes very difficult because of the termi- 
nating Gaussian hypergeometric series 2F1 in ()6.12p (prob- 
ably, this would be a nontrivial research project in its own 
right). However, we can argue that the function exp(uz) 
with /z, M, 2: G C is only analytic at z = in the sense of 
complex analysis if is a nonnegative integer, fi = m with 
m G No, yielding z"'' exp{uz) = u"z™+"/n!. If jj. is 

nonintegral, yU G C \ No, a power series around z — does 
not exist. 

Thus, we can conclude that also for /? / 7, the one- 
center limit r' = of the rearranged addition theorems 
(|6.8p for XN.hiP^''' ^ '"'') "ioss not exist if the principal 
quantum number N is nonintegral, G M \ N. 

From a mathematical point of view, a one-range ad- 
dition theorem for a function /(r ± r') is a mapping 
R-^ X C. We are interested in one-range addition 

theorems because we would like to have a unique (infinite 
series) representation of f{r±r') with separated variables 
r and r' that is valid for the whole argument set M'^ x M''. If 
we accept this premise, then we have to conclude that Gu- 
seinov's manipulations, which produced the rearranged ad- 
dition theorems (j6.8p for Xn Lil^^'''^'''')^ ^^'^ least in the 
case of nonintegral principal quantum numbers A'^ G M \ N 
a complete failure. 

This observation does not rule out the possibility that 
an appropriate reinterpretation might make the rearranged 
addition theorems ()6.8p with G M \ N mathematically 
meaningful in a restricted sense as an approximation, al- 
though it does exist for the whole argument set M'^ x R^. 
This has to be investigated. Moreover, it is still open 
whether the rearranged addition theorems (|6.8p exists for 
the whole argument set R^ x R'^ if the principal quantum 
number N of the Slater-type function is a positive integer, 

G N. Also this remains has to be investigated. In all 
cases, the burden of proof lies with Guseinov. 

These observations should be particularly worrisome 
for those who advocate the use of Slater-type functions 
with nonintegral principal quantum numbers as basis func- 
tions in molecular electronic structure calculations. It is 



generally accepted that a good basis set should produce 
highly accurate approximations, but it should also lead to 
multicenter integrals that can be computed efficiently and 
accurately at tolerable computational costs. These two 
requirements have so far been mutually exclusive. 

Slater-type functions with integral principal quantum 
numbers or other exponentially decaying function sets are 
well suited to produce highly accurate approximations 
in electronic structure calculations, but their multicen- 
ter integrals are notoriously difficult. In contrast, Gaus- 
sian functions are nonphysical. Consequently, large ba- 
sis sets are needed to accomplish sufficiently accurate re- 
sults. However, Gaussian functions have one, albeit deci- 
sive advantage: Their multicenter integrals can be com- 
puted comparatively easily. 

If we nevertheless want to use physically better mo- 
tivated exponentially decaying basis functions, we should 
concentrate on those functions that promise the simplest 
multicenter integrals, even if we should have to sacrifice 
some accuracy. Superficially, Slater-type functions with 
nonintegral principal quantum numbers look attractive 
since they promise somewhat better results than Slater- 
type functions with integral principal quantum numbers. 
But they achieve this at the cost of significantly more 
complicated multicenter integrals. Therefore, we would 
first need some fundamental mathematical breakthroughs 
to make Slater-type functions with nonintegral principal 
quantum numbers competitive with Slater-type functions 
with integral principal quantum numbers or with other ex- 
ponentially decaying functions. As far as I can judge it, 
these breakthroughs are not in sight. 

There is another, more general remark which I would 
like to make. Addition theorems of exponentially decay- 
ing functions as well as their multicenter integrals are fairly 
complicated mathematical objects. Accordingly, it is usu- 
ally very difficult and often even practically impossible to 
justify nontrivial manipulations by rigorous mathematical 
proofs. This is highly embarrassing. The only way out, 
which I see, is to be very conservative and cautious in or- 
der to be on the safe side. 

Let us assume that / belongs to a suitable Hilbert 
space. If we expand f{r ± r') in terms of a function set, 
which is complete and orthogonal in that Hilbert space, 
then the resulting one-range addition theorem converges 
in the mean with respect to the norm of the correspond- 
ing Hilbert space. Moreover, in many cases the Schwarz 
inequality (see for example 82, Eq. (6) on p. 31]) suffices 
to guarantee that the use of this addition theorem in a 
multicenter integral produces a convergent expansion. 

However, we are no longer on the safe side if we re- 
place the complete and orthogonal functions by other func- 
tions that are only complete but not orthogonal. For- 
mally, this may yield another one-range addition theorem 
for /(r ± r'), but we do not know whether this expansion 
makes sense or not. This has to be proved. If we cannot 
prove this, then we should rather avoid such a possibly 
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dangerous manipulation. It is not acceptable to ignore 
potential problems of that kind and hope for the best, as 
it was done by Guscinov. 

7 Weakly Convergent One-Range 
Addition Theorems for the 
Coulomb Potential 

Guseinov derived one-range addition theorems for the 
Coulomb potential via addition theorems for the Yukawa 
potential exp(— /3r)/r [76j . which can be viewed to be an 
exponentially screened Coulomb potential and which is 
also a special Slater-type function according to 

^ = (47r)i/2/3x%(/3,r). (7.1) 

Guseinov's idea was to set iV = L = M = in his un- 
symmetrical one-range addition theorems for Slater-type 
functions discussed in Section [S] and to perform the limit 
/3 — > 0. Unfortunately, things are more complicated and 
Guseinov's simplistic approach is fundamentally flawed. 

For fc = 0, 1, 2, . . . , the Yukawa potential belongs to the 
weighted Hilbert space Ll^(M.^) defined by but not 

for k = —1 (see for example [ini p. 410]). This was appar- 
ently overlooked by Guseinov who used his functions with 
unspecified k. Accordingly, the addition theorems (|6.1|1 
with N — L — M = converge for fc = 0, 1, 2, . . . in the 
mean with respect to the norm (|4.19p of ^^^(R^), but not 
for fc = — 1. 

Guseinov replaced in his one-range addition theorems 
(|6.ip the overlap integrals (j6.1b[) involving his complete 
and orthonormal functions k'^^eil,''') by overlap inte- 
grals of Slater-type functions according to ()6.5p . This 
step, which yields the addition theorems ()6.6p . is math- 
ematically legitimate, but numerically dubious. Also the 
next step - the replacement of the expansion functions 
k^^gi'y,'!') on the right-hand side of (|6.6p by Slater-type 
functions according to (|6.4p , which yields the addition the- 
orems ()6.7p - is in the case of the Yukawa potential at least 
for fc = 0, 1, 2 . . . mathematically justified, but only as long 
as the order of the infinite summations in (j6.7p is retained. 

Guseinov's final step - the rearrangement of the infinite 
summations in (|6.7p which formally produces the addition 
theorems ()6.8p with Slater-type functions as expansion 
functions - is not legitimate: The one-center limit of these 
addition theorems does not exist for N = L = M = Q: 
If we set r' = in the addition theorems ()6.ip with 
N = L = M = 0, we obtain expansions that are special 
cases of the Laguerre expansion (|6.12p for exp(uz) with 
fj, = —I. This Laguerre expansion cannot be reformulated 
as a power series in z since exp(wz)/z is not analytic at 
z ~ in the sense of complex analysis. Accordingly, the 
rearranged addition theorem (|6.8p with N = L = M = 
does not exist for the whole argument set x R'^. 



Thus, Guseinov's derivation [39^ of one-range addition 
theorems for the Coulomb potential with Slater-type func- 
tions as expansion functions fails since his starting point 
- the unsymmetrical one-range addition theorems (|6.8p 
with A^ = L = M = 0-do not exist for all arguments 
r,r' e R3. 

As a possible remedy, we could expand the Coulomb 
potential in terms of a suitable complete and orthonor- 
mal set of functions. If we expand l/\r — r'\ in terms 
of Guseinov's functions, we formally obtain the following 
symmetrical one-range addition theorems which are spe- 
cial cases of the general addition theorem (|5.22p : 

nim 
n i'm' 

X ^^^/P,r) ,^^,[AP,r'), (7.2a) 

n'i'm'xP) 

= I (/3,0]*^'%C:::,(/3,r')dV, (7.2b) 

,C-,(/3,r') = J (/3,r-)]*^^dV. (7.2c) 

Analogous symmetrical addition theorems can also be de- 
rived for the Yukawa potential. 

Instead of Guseinov's functions, we could use in (|7.2p 
or in analogous addition theorems for the Yukawa po- 
tential any other function set that is complete and or- 
thonormal in an appropriate Hilbert space. An obvious, 
albeit Gaussian-type choice would be the eigenfunctions 
r2™^(/3,r) of the three-dimensional isotropic harmonic os- 
cillator, which are defined by (j4.7p and which are complete 
and orthonormal in L^(R'^). 

In the case of ()7.2p there is a principal problem that was 
apparently overlooked by Guseinov: The Coulomb poten- 
tial does not belong to any of the Hilbert spaces considered 
in this article or implicitly used by Guseinov, since they 
all involve an integration over the whole R'^. Thus, for 
k = —1,0,1,2,... the addition theorems (j7.2p diverge in 
the mean with respect to the norm ()4.19p of L^f, (R'^). 

It is tempting to conclude that all attempts of con- 
structing one-range addition theorems for l/|r — r'| by ex- 
panding it in terms of function sets, that are complete and 
orthonormal with respect to an inner product involving 
an integration over the whole R^, are doomed. However, 
this conclusion is premature, and expansions of the type 
of (|7.2p may well be our best chance of achieving our aim. 

As discussed in Section [51 one-range addition theorems 
are based on Hilbert space theory and thus rely heavily on 
concepts from approximation theory. Loosely speaking, 
we may say that it is the purpose of approximation theory 
to provide convenient expressions - for example series ex- 
pansions " that allow an efficient and reliable evaluation 
of a given mathematical object. 

There is, however, a very important difference between 
the use of one-range addition theorems in multicenter in- 
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tegrals and more conventional applications of approxima- 
tion theory: We are not interested in evaluating the func- 
tions that are represented by addition theorems. We only 
use these addition theorems to simplify the integrations in 
multicenter integral. Addition theorems are only interme- 
diate results which ultimately produce series expansions 
for multicenter integrals. 

The use of convergent expansions in integrals has the 
undeniable advantage that normally only comparatively 
mild assumptions are needed to guarantee that integration 
and summation can be interchanged and that the resulting 
expansions converge. Nevertheless, the use of convergent 
expansions in integrals is to some extent a luxury and not 
strictly necessary. We are free to use a divergent expansion 
in an integral and interchange integration and summation 
if we can guarantee that the resulting expansion converges 
to the correct result. 

Obviously, such an approach gives us additional pos- 
sibilities, but it would be naive to expect a free lunch: 
It is grossly negligent to use divergent series in integrals 
without knowing criteria of manageable complexity that 
guarantee the convergence of the resulting expansions. 

As discussed in Section [l] two-range addition theo- 
rems of the type of the Laplace expansion ()1.2p con- 
verge pointwise since they are essentially rearranged three- 
dimensional Taylor series in the Cartesian components of 
the shift vector. By relaxing the requirement of point- 
wise convergence, we arrive at one-range addition theo- 
rems which in general converge in the mean or strongly 
with respect to the norm of a suitable Hilbert space. 

As a consequence, one-range addition theorems for 
functions, which do not belong to a given Hilbert space, 
diverge with respect to its corresponding norm. If we nev- 
ertheless insist on extending the formalism of one-range 
addition theorems to functions that do not belong to this 
Hilbert space, the requirement of convergence in the mean 
or strong convergence is too demanding. We have to re- 
place it by a weaker type of converge that is related to 
and inspired by the theory of generalized functions in the 
sense of Schwartz [129] . 

Quantum theory is largely based on functionals, which 
are mappings from a function space to the complex num- 
bers. Accordingly, the numerical values and the local prop- 
erties of wave functions are - strictly speaking - irrelevant, 
and convergence of functionals, also called weak conver- 
gence, is all we really need. Stronger types of convergence 
such as pointwise convergence or convergence in the mean 
are in principle pure luxury. 

Thus, the one-range addition theorems (|7.2p . which di- 
verge for any fc = — 1,0,1,2,... with respect to the norm 
(|4.19p of the Hilbert space L^fc(]R^), should be interpreted 
as weakly convergent expansions that produces correct re- 
sults when used in suitable functionals. 

In the context of multicenter integrals, weakly conver- 
gent expansions have largely been a terra incognita, but 
there is a detailed mathematical literature both on weak 



convergence in general (see for example [130' Chapters 10 
and 11]) as well as on the expansion of generalized func- 
tions or distributions in terms of orthogonal polynomials 
(see for example [1311 11321 1133] and references therein). 

Weakly convergent expansions can be quite useful. In 
|18| weakly convergent expansions for exp(ip • r) in terms 
of complete orthonormal and biorthogonal function sets 
were constructed. In some cases, these expansions simplify 
the evaluation of Fourier transforms, and they can also be 
used for the construction of one-range addition theorems 
(see m Section VII] or Q^). 

Expansions for the plane wave, that closely resemble 
those derived in [TH], were also constructed by Guseinov 
i28j Eqs. (45) - (46)]. Guseinov, who did not mention the 
article [IH] in [5S] , either overlooked or ignored the obvious 
fact that that the plane wave does not belong to any of 
the Hilbert spaces which he implicitly used. Accordingly, 
Guseinov's expansion diverge in the mean and can only 
converge weakly. Guseinov's oversight is hard to under- 
stand because he had cited [Hj in several other articles 
|25l[Ml[3l[40l|4ll|49l[52l[53l[57]. 

In my opinion, the essential features of weak conver- 
gence in contrast to strong convergence can be explained 
most easily via the Euclidean vector space C°° of infi- 
nite row or column vectors u — (ui,U2, ■ ■ ■) with com- 
plex coefficients m„. By equipping the vector space C°° 
with the inner product {u\v) = X^^^i I^"]*^"; o^)- 
tain the corresponding Hilbert space C C°° with norm 
]]it]| = ^/(u\u). The condition \\u\\ < oo, which defines 
the Hilbert space £'^, can only be satisfied if the coefficients 
Un € C of tt decay sufficiently fast as n — s- oo. 

Let us now assume that some vector w e C°° cannot 
be normahzed, i.e., \\w\\^ = {w\w) = J2^=i l^nj^ < °° 
does not hold. Thus, w ^ £'^, but this does not imply that 
all inner products {w,u) with normalizable u € £^ do not 
exist. If the coefficients m„ of u decay sufficiently rapidly 
as n ^ cxo, then the infinite series (w\u) = J2'^=ib^n]*Un 
may well converge to a finite result, although the series 
{w\w) ^ V^ri? diverges. 

It is important to note that |(i(;|m)| < oo cannot hold 
for arbitrary u. G but only for a suitably restricted sub- 
space yy C C of normalizable vectors such that 
a; G W implies |(i(;|a;)| < oo. Accordingly, my interpreta- 
tion of weak convergence resembles at least conceptually 
the theory of rigged Hilbert spaces or Gelfand triplets (see 
for example |1341 I135j and references therein). A very 
readable account of rigged Hilbert spaces from the per- 
spective of quantum mechanics and their relationship with 
Dirac's bra and ket formalism can be found in the book 
by Ballentine [136 . Chapter 1.4]. 

In principle, it is desirable to specify for a given w the 
whole set W, but this may be difficult. In practice, it may 
well be sufficient to specify only a sufficiently large subset 
of W. Let us for instance assume that the coefficients Wn 
of a vector w ^ £j^ are all finite and satisfy ]i(j„] ^ with 
> —1/2 as n — » oo. Thus, w cannot be normalized. 
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If, however, the coefficients w„ of a vector a; S £^ are all 
finite and satisfy \LUn\ ^ n" with a < —(3 — 1 as n — > oo, 
then the infinite series (tclw) = Yl'^=i{^'n\*^n converges 
and the inner product (if |<^) makes sense. 

However, the condition |w„| n" with a < —[3 — I 
as n — * cx) does not suffice to specify the whole set W of 
vectors us with |(u;|a;)| < oo. It can happen that the se- 
ries Yl^=i['^ri]*^n diverges but is summable. If summabil- 
ity techniques are included in our arsenal of mathematical 
techniques, we could even discard the otherwise essential 
requirement a; G £^ and try to make inner products {w\'jj) 
with both w,u) ^ mathematically meaningful. A (very) 
condensed review of the classic summability methods as- 
sociated with the names of Cesaro, Abel, and Riesz can be 
found in Zayed's book [1371 Chapter 1.11.1]. More detailed 
treatments of linear summability methods can be found in 
specialized monographs such as the books by Boos [ISSj . 
Hardy [HS], Knopp [HOI, and Powell and Shah [HT] . 

An interpretation of one-range addition theorems of 
the type of (|7.2p as weakly convergent expansions requires 
that we first specify the functional, in which addition the- 
orems are to be used. Then, regularity criteria have to 
be formulated whose validity guarantee the convergence 
of the resulting expansions for this functional. 

The most difficult integrals, which occur in the LCAO- 
MO approach, are the notorious two-electron integrals 
C{f,g) defined by (|6.3p . whose evaluation can become ex- 
tremely difficult if / and g are nonclassical two-center 
charge distributions of exponentially decaying functions. 

A theoretical analysis of the integrals C{f,g) is by no 
means easy (see for example [1421 Chapter 9]). In this 
context, it is instructive to replace in (|6.3p the Coulomb 
potential by the Yukawa potential, yielding 

W.,,«.///>)£^^Z^a=.dV ,7.3, 

Obviously, we have C{f,g) = Hm^^o 3^(/j 5; /?)■ Unfortu- 
nately, it is not guaranteed that this limiting process is 
continuous and produces a finite result. It can be shown 
with the help of Fourier transformation (see for example 
|1431 Section II] and references therein) that 3^(/, g\ (3) ex- 
ists for arbitrary densities /, 5 G L^(K^), but this does not 
suffice to guarantee the existence of C{f,g). The densi- 
ties / and g have to satisfy more sophisticated conditions 
than square integrability to guarantee this (see for example 
|1421 Section 4.3 (Hardy-Littlewood-Sobolev inequality)] 
or |144| Example 3 (Sobolev's inequality) on p. 31]). 

If we want to apply weakly convergent addition the- 
orems of the type of (|7.2p for the evaluation of two- 
electron integrals, we need additional criteria that guar- 
antee the convergence of the resulting expansions to the 
value of C{f,g) for physically and mathematically reason- 
able charge densities / and g. 

If we insert the addition theorems (|7.2|) into the inte- 



gral (|6.3p . we obtain the following series expansions: 

= Y. .r:;f?'„,(ffljF:;,(s)tG:;;,,(/3). (7,4.,) 

nim 
n i'ni 

feF™,(/3) = J r{r) ,^^/p,r)d^r, (7.4b) 

,G;V(/3) - !<?(/) feVl/;";,,(/3,r')dV. (7.4c) 

If only a finite number of the integrals j,F™^(/3) and 

j.G™£,(/3) are nonzero, then the use of (|7.2p in C{f,g) 
obviously produces the correct result. 

In general, there will be infinitely many nonzero inte- 
grals fcF™^(/3) and j,G™^,(/3). In this case, asymptotic 
estimates for the coefficients j,r"f™^, (/?) defined by ()7.2bp 
have to be constructed that hold in the case of large in- 
dices. Subsequently, asymptotic conditions on the decay of 
the integrals fcr™^(/?) and ^G™, £/(/?), respectively, in the 
case of large indices can be formulated which guarantee 
the convergence of the series expansion ()7.4p . 

It cannot be denied that such a proof of the conver- 
gence of the series expansions (17. 4p - although most likely 
manageable - would be highly pedestrian. Nontrivial tech- 
nical difficulties are also quite likely. Obviously, a more 
elegant convergence proof based on concepts of functional 
analysis and suitable functions spaces would be highly de- 
sirable. Unfortunately, this is currently not in sight. 

There is also another problems. Kato 145J showed that 
the singularities of the potential of atomic and molecular 
Hamiltonians produce discontinuities of the wave functions 
commonly called Coulomb or correlation cusps (see also 
jl46 ]). This implies that atomic or molecular wave func- 
tions do not possess continuous derivatives of arbitrary or- 
der at the locations of the nuclei. Consequently, exponen- 
tially decaying functions, which cannot be differentiated 
arbitrarily often at the origin, provide much better ap- 
proximations than infinitely differentiable Gaussian-type 
functions. According to experience, a few exponentially 
decaying functions normally suffice to accurately model 
the discontinuities of atomic and molecular wave functions 
at the nuclei. But there is a price: Exponentially decay- 
ing functions are fairly complicated mathematical objects 
whose multicenter integrals are notoriously difficult. 

It is quite likely that we also have to pay some price 
if we try to analyze the weak convergence of the addi- 
tion theorems ()7.2p in appropriate functionals such as the 
integral ()6.3p . As remarked above, weak convergence is re- 
lated to Schwartz's theory of generalized functions which 
applies to integrals containing a distribution multiplied by 
a so-called test function. To make these integrals well de- 
fined, the test functions must possess a sufficient amount 
of mathematical "niceness" in order to compensate the 
possibly highly irregular behavior of generalized functions 
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such as the Dirac delta "function" or their derivatives. 

In mathematics, it is common to use test functions 
with highly idealized properties because this greatly facil- 
itates proofs. A frequently used test function space is the 
Schwartz space iS(M'^) of rapidly decreasing and infinitely 
differentiable functions ip : M."^ C, which satisfy 



sup 



d_ 

dx 



< oo 



(7.5) 

for all k,m,n,u,v,w € No (see for example [HI p. 133]). 
The oscillator eigenfunctions il™^(/3, r) defined by (j4.7p or 
also other Gaussian-type functions belong to 5(M'^), but 
none of the exponentially decaying function sets consid- 
ered in this article. 

Even more idealized are test function tp: M.^ C be- 
longing to the space I?(M'^) consisting of infinitely differ- 
entiable functions with compact support (see for example 
dm Chapter 6.2]). 

Because of their infinite differentiability, test func- 
tions belonging to 5(R'^) or I?(IR'^), respectively, are ex- 
tremely convenient analytical tools. Their use greatly sim- 
ply mathematical proofs. However, infinitely differentiable 
functions are not ideally suited for the representation of ef- 
fective one-particle wave functions which are discontinuous 
and also decay exponentially (see for example [147|, I148j 
and references therein). 

The completeness of the oscillator eigenfunctions 
il™^(/3,r) in L^(E'^) implies that they can approximate 
effective one-particle wave functions in the mean with an 
in principle unlimited precision, but one should not expect 
rapid convergence. If we instead use physically better mo- 
tivated exponentially decaying functions, which are not 
infinitely differentiable, we can expect faster convergence, 
but all mathematical manipulations involving these func- 
tions will be (much) more difficult. Most likely, this also 
applies to the analysis of the properties of weakly conver- 
gent addition theorems of the type of (|7.2p . 

Both the index k as well as the scaling parameter (3 of 
the expansion functions fc^^f(/?, J") influence the rate of 
convergence of the infinite series in ()7.4ap . It is an obvious 
idea to try to maximize convergence by optimizing both 
k and p. Unfortunately, our current level of understand- 
ing does not allow detailed predictions. Most likely, the 
"optimal" values of k and /3 depend quite strongly on the 
properties of the two charge densities / and g. It is also 
possible that for given / and g only some values of k will 
lead to a convergent series expansion (|7.4ap . 

The use of a nonzero k can at least potentially lead to 
stabihty problems. For fc = 0, the integrals feF™^(/3) and 

j,G™ £/(/?) are the coefficients of expansions of / and g* in 
terms of Lambda functions, which are according to (j4.10p 
orthonormal with respect to an integration over the whole 
R"^. Since the coefficients of orthogonal expansions sat- 
isfy Parseval's equality (|2.5p . the integrals j,F™^(/3) and 

j.G™^/(/3) with fc = are bounded in magnitude and 



vanish for large indices. For fc 7^ 0, Guseinov's functions 
fe^™£(/^' r) are not orthogonal with respect to an integra- 
tion over the whole M^, but satisfy ([iTfll . Thus, fcF;^^(/3) 

and fcG™ ^/(/3) with fc 7^ are not the coefficients of or- 
thogonal expansions and do not satisfy Parseval's equality. 
Accordingly, it is not guaranteed that they are bounded in 
magnitude and that they vanish for large indices. 

Let me emphasize that the convergence or divergence 
of a one-range addition theorem and the convergence or 
divergence of the resulting expansion for a multicenter in- 
tegral are not directly related. In [39l pp. 212 - 213] Gu- 
seinov claimed that the convergence of his addition theo- 
rems for the Coulomb potential can be demonstrated via 
the convergence of the resulting series expansions for three- 
center integrals of the type of (|6.2p . This is of course 
wrong: Weak convergence does not imply convergence in 
the mean, let alone pointwise convergence (see for instance 
[149, §29]). Moreover, Guseinov should know that an ob- 
served agreement of different floating point computations 
up to a certain number of digits does not necessarily prove 
anything (see for example |150, Table 2]). 

But even if the one-range addition theorems (|7.2p or 
related addition theorems should lead to divergent series 
expansions for multicenter integrals, the resulting expan- 
sions can nevertheless be computationally useful. In the 
book by Bornemann, Laurie, Wagon, and Waldvogel [1511 
p. 225], there is the following instructive remark: 

The question whether a series converges is 
largely irrelevant when the reason for using a 
series is to approximate its sum numerically. 

Thus, even wildly divergent series can be used for compu- 
tational purposes if suitable summation methods can be 
found. 

I am aware of several predominantly theoretical articles 
dealing with the summation of divergent Laguerre expan- 
sions (see for example |152[ 1153") and references therein). 
However, from a practical point of view, it is probably 
more effective to use purely numerical summation meth- 
ods. I am also skeptical about the classical summabil- 
ity methods mentioned above and discussed in books by 
Boos [138], Hardy [l39], Knopp [ITQ;, and Powell and Shah 
j l41| . Instead, I suggest to use nonlinear sequence trans- 
formations, which can also be applied in the case of slowly 
convergent sequences and series and which often achieve 
spectacular improvements of convergence. In that con- 
text, it may be of interest to note that the most recent 
(third) edition of the book Numerical Recipes 154] now 
also discusses nonlinear sequence transformations (for fur- 
ther details, see also |155j ). 

The best known and most often used sequence transfor- 
mations are the so-called Pade approximants |156| which 
accomplish an acceleration of convergence or a summation 
by converting the partial sums of a power series to a dou- 
bly indexed sequence of rational functions. As documented 
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by the long list of successful applications in the book by 
Baker and Graves-Morris 157J, Pade approximants have 
become the standard tool in theoretical physics and in ap- 
plied mathematics to overcome problems with slowly con- 
vergent or divergent power series. 

It is, however, not so well known among non-specialists 
that alternative sequence transformations as for example 
Wynn's epsilon and rho algorithm 158^ 1159] . Brezinski's 
theta algorithm [160j . or Levin's transformation |161j and 
later generalizations (see for example [1621 11631 1164] and 
references therein) can at least for certain computational 
problems be much more effective than Pade approximants. 

As emphasized already several times, we have no a pri- 
ori reason to assume that one-range addition theorems 
necessarily lead to rapidly convergent expansions for mul- 
ticenter integrals. Accordingly, our ability of evaluating 
infinite series effectively and reliably is crucial for the prac- 
tical usefulness of addition theorems. 

The conventional process of adding up the terms of a 
series successively is at least in principle able to produce 
approximations with unlimited accuracy for a convergent 
series, but it completely fails in the case of a divergent 
series. Moreover, adding up the terms successively is in 
far too many cases prohibitively inefficient. Therefore, it 
makes sense to try to use sequence transformations for 
the evaluation of series expansions for multicenter integrals 
whenever possible. This is not a new idea. The oldest ar- 
ticle using sequence transformations for the evaluation of 
multicenter integrals, which I am aware of, was published 
in 1967 by Petersson and McKoy [165j . 

Personally, I became interested in sequence transforma- 
tions during my PhD thesis [120] . in which series expan- 
sions for multicenter integrals played a major role. Since it 
(too) often happened that my series expansions converged 
slowly (see for example |1191 Table II]), it was a natural 
idea to speed up convergence with the help of sequence 
transformations. During my PhD thesis, I only knew lin- 
ear series transformations as described in the classic, but 
now outdated book by Knopp 140J, which turned out to 
be ineffective. I was completely ignorant of the more pow- 
erful nonlinear transformations, which often accomplish 
spectacular improvements of convergence. Only later, I 
used nonlinear sequence transformations with considerable 
success for the evaluation of multicenter integrals of expo- 
nentially decaying functions [IMl [Ml [Ml [Ml [HO] . 

The usefulness of sequence transformations in the con- 
text of multicenter integrals is not limited to the evaluation 
of infinite series representations. Particularly noteworthy 
seems to be Safouhi's approach. Starting from the Fourier 
transform (|5.3p of B functions, Safouhi converts compli- 
cated multicenter integrals of B or Slater-type functions 
to multidimensional integral representations that have to 
be evaluated by numerical quadrature. 

At first sight, this does not look like a good idea be- 
cause the oscillatory nature of the multidimensional inte- 
gral representations makes the straightforward application 



of conventional quadrature methods difficult. However, 
these problems can be overcome by combining quadrature 
schemes with suitable nonlinear sequence transformations. 
Based on previous work of Sidi 1171' and of Levin and Sidi 
[.172] , Safouhi succeeded in developing some extrapolation 
techniques specially suited to his needs. This permits a re- 
markably efficient and reliable evaluation of complicated 
multicenter integrals via oscillatory integral representa- 
tions (see for example [T73 [HI [Ml [Hi [113 [Hi fT75] 
and references therein). 

In my opinion, Safouhi's work is a convincing demon- 
stration of the practical usefulness of extrapolation and 
convergence acceleration techniques in electronic structure 
calculations. Safouhi's work also shows that there is no 
reason for despair if existing sequence transformations turn 
out to be not powerful enough to solve certain problems. 
It may well be possible to construct new transformations 
that can do the job. 

Sequence transformations can be useful also in com- 
pletely different contexts. For example, I have applied se- 
quence transformations successfully in such diverse fields 
as the evaluation of special functions and related objects 

[Ml [Ml [HH IMl [Ml IMl [Ml IMl [Ml [Ml [M], 

the summation of strongly divergent quantum mechani- 
cal perturbation expansions [TM', IW, IHI [Ml [T^. TMl 
,195. 182> ITM [TM [Ml [Ml il99, 2001 [Ml [^ 203j, 
the prediction of unknown series coefficients [1901 11941 
11951 12011 1204j , and the extrapolation of quantum chemical 
crystal orbital and cluster electronic structure calculations 
for oligomers to their infinite chain limits of stereoregu- 
lar gwasi-onedimensional organic polymers [2051 12061 1207] . 
Many other applications of sequence transformations are 
listed in my Habilitation thesis |184j and in |164j . 

These examples should convince even a skeptical reader 
that sequence transformations are extremely useful com- 
putational tools, and that it is worth while to invest 
time and effort to understand their power as well as their 
shortcomings. Of course, sequence transformations are no 
panacea, and not all my attempts were successful. How- 
ever, even my failures often turned out to be fruitful in 
the long run, since they frequently provided new insight 
which ultimately paved the way for the construction of 
new convergence acceleration and summation techniques 
[163, 164, 167, 183, 186, 187l [T96l[208l[209| . 

In my opinion, it is hard to understand that there are 
still researchers like Guseinov who try to evaluate series 
expansions for multicenter integrals without the help of se- 
quence transformations. Multicenter integrals of exponen- 
tially decaying functions are notoriously complicated, and 
a nice explicit expression for a multicenter integral does 
not necessarily permit its efficient and reliable evaluation. 
It is also necessary to employ powerful and sophisticated 
numerical techniques. The conventional process of adding 
up the terms of a series successively until convergence is 
achieved does not fall into this category. 

The currently most complete reference on Pade ap- 
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proximants is the impressive book by Baker and Graves- 
Morris [157] . The most recent monograph on sequence 
transformations is the book by Sidi [210] . It contains a 
wealth of information and is undeniably useful for special- 
ists, but presentation and choice of topics makes it basi- 
cally unsuited for novices (compare also my book reviews 
|2111 1212] ). So, the best choice may well be the book by 
Brezinski and Redivo Zaglia |213j . Some people also like 
my long review |163j . As a first introduction for novices, I 
can recommend Appendix A of the book by Borneniann, 
Laurie, Wagon, and Waldvogel |,151j . 

8 Differentiation Techniques and 
the Spherical Tensor Gradient 
Operator 

Often, it comparatively easy to obtain explicit analytical 
expressions for multicenter integrals over scalar functions, 
which are irreducible spherical tensors of rank zero with 
respect to their local (atomic) coordinate systems. If, how- 
ever, the functions occurring in the multicenter integral are 
irreducible spherical tensors of higher ranks, the derivation 
of explicit expressions can easily become extremely diffi- 
cult or even practically impossible. 

It is also an empirical fact that it is usually much eas- 
ier to differentiate than to integrate (integration is an art, 
but every fool and in particular also computer algebra sys- 
tems like Maple or Mathematica can differentiate). Ac- 
cordingly, it is an obvious idea to try to generate an ex- 
plicit expression for a difficult multicenter integral over 
nonscalar functions by differentiating the simpler multi- 
center integral over scalar functions - preferably the sim- 
plest scalar functions - with respect to scaling parameters 
and/or nuclear coordinates (see for example |12H Section 
IV]). This approach is not restricted to the notoriously dif- 
ficult multicenter integrals of exponentially decaying func- 
tions. It is also used in the case of multicenter integrals of 
the so-called spherical Gaussian functions (see for example 
[aaiHIlEIllElTlEIl] and references therein). 

It is relatively easy to generate multicenter integrals 
of higher scalar functions by differentiating the simplest 
scalar functions with respect to their scaling parameters. 
In the case of a Is Slater-type or Gaussian function, we 
can construct higher scalar functions easily by repeatedly 
using the relationships d exp{—ar)/da = — rexp(— ar) or 
d exp{—ar^) / da = — exp(— ar^). 

The generation of anisotropic functions, which are ir- 
reducible spherical tensors of rank £, from scalar functions 
is less straightforward, but also here an highly developed 
mathematical technology based on the differential opera- 
tor 3^™(^) available (see for example 3 and [4, Section 
3] and references therein). 

In Section 3], Guseinov used differentiation tech- 
niques for the derivation of more complicated addition 



theorems by differentiating his one-range addition theo- 
rems for l/\r — r'\ with respect to the Cartesian compo- 
nents of either r — {x,y,z) or r' = (x',y',z'). Guseinov 
pursued essentially analogous approaches in his articles 

[3iiEa[3B,^^i38i^42,43,45.d6i[saiiaEniiii]. 

It is not a new idea to generate more complicated ad- 
dition theorems of anisotropic functions by differentiating 
simpler addition theorems of scalar functions. In ^219,, Sec- 
tions IV and VI] , it was shown that the two-range addition 
theorems of the irregular solid harmonic and the modified 
Helmholtz harmonic, respectively, can be derived by differ- 
entiating the two-range addition theorems for the Coulomb 
potential and the Yukawa potential, respectively. Similar 
ideas were pursued in ll6 . 

Unfortunately, this does not imply that Guseinov's ap- 
proach is computationally efficient or at least mathemat- 
ically sound. This Section first describes how irreducible 
spherical tensors can be differentiated comparatively eas- 
ily with respect to the Cartesian components of their ar- 
gument vectors. The second topic is less technical and 
refers to more fundamental mathematical problems: In 

[HI 132 [S3 Ell [33 [Ml [Ml HIl ia 113 HSl ia IMl El EHl ^ 

Guseinov consistently ignored all questions of convergence 
and existence of his addition theorems obtained by differ- 
entiating simpler addition theorems. 

This article considers exclusively addition theorems of 
irreducible spherical tensors of the type of (jl.ip . More- 
over, all functions occurring in these addition theorems 
are also irreducible spherical tensors of a given rank ^ > 0. 
Irreducible spherical tensors F™(r ± r'), that are of in- 
terest in the context of electronic structure calculations, 
can all be differentiated at least a finite number of times 
with respect to the Cartesian components of either r or 
r' . Because of the convenient orthonormality properties 
of the spherical harmonics, it is highly desirable to express 
the angular part of Cartesian differentiations in terms of 
spherical harmonics. The necessary algebra can be done, 
but in particular for large angular momentum quantum 
numbers of the irreducible spherical tensors and for large 
orders of the differential operators, we would be confronted 
with messy expressions and nontrivial technical problems. 
There is also the danger that delta function contributions, 
which can occur in multicenter integrals and addition the- 
orems (see for example the articles by Pitzer Kern, and 
Liscomb 220 and by Kay, Todd, and Silverstone [221j or 
the book by Judd [222', Chapter 5.3]), are easily overlooked 
if we differentiate with respect to Cartesian components 
(see also PSII^I^I^ V 

A much more convenient approach is possible that com- 
pletely avoids all differentiations of irreducible spherical 
tensors f ™(r ± r') with respect to Cartesian components 
of the argument vectors and only requires differentiations 
with respect to the radial variables r = |r| and r' = \r'\, 
respectively. This is accomplished by reformulating the 
original differential operators with unspecified transforma- 
tion properties as finite sums of differential operators that 
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are irreducible spherical tensors of a given integral rank. 

Let us assume that Vn{r) is a polynomial of degree n 
in the Cartesian components of r = {x,y, z): 



u-\-v^w<n 



z'^ 

uvw y 



(8.1) 



The regular solid harmonic 3^™('') is according to (|A.4p a 
homogeneous polynomial of degree I in the Cartesian com- 
ponents of r = (x,y,z). Thus, the completeness and the 
orthonormality of the surface spherical harmonic Y^{9, 0) 
with respect to an integration over the unit sphere in 
implies that the polynomial Vnir) can be expressed as a fi- 
nite sum of solid harmonics in r multiplied by even powers 
oir=\r\ (see also EUl §96 on pp. 147 - 148]): 



2i/+A<n 

E 



X 

E 



y'xir). 



(8.2) 



This is a relationship among polynomials in the Cartesian 
components of an essentially arbitrary three-dimensional 
vector r. Thus, ()8.2p also holds if we replace r by V = 
{d/dx,d/dy,d/dz), and any differential operator, which 
is polynomial of degree n in the Cartesian components of 
V, can be expressed as a finite sum of products of integral 
powers of the Laplacian = d'^/dx'^ + d'^/dy'^ + d'^/dz'^ 
multiplied by a solid harmonic in V: 



U+V+W<n / ^ \ U / \ V / n 

V c(") —] -] - 

u.v,w>() ^ ^ \ y/ \ 



2u+X<n X 

= E E Cl"iv-3^nv). 

i^,X>0 fJ.=-X 



(8.3) 



We obtain an explicit expression for the spherical ten- 
sor gradient operator y™(V) by replacing in (IA.4P the 
Cartesian components of r by those of V: 



2£+l 
An 



{i + my.{£ - my. 



1/2 



\ dx dy J \dx dy J \ dz ^ 

2™+2fe(TO + k)\k\{£ - m - 2fc)! 



E 



(8.4) 



In Martin's book [4, p. 62], yp{V) is called the Erdelyt 
operator. I prefer instead the name spherical tensor 
gradient operator because the expressions for products 
y^^'' {'^)Fl^^ {r), that will be mentioned later, generalize - 
as emphasized by Dayman |228| p. 2558] - the well known 
gradient formula in angular momentum theory (see for ex- 
ample [m Chapters 5.7 and 5.9] or [230, Chapter 11.11]). 

In Martin's book [H Section 3] , the mathematical the- 
ory of y^^iV) as well as numerous applications predomi- 
nantly in classical physics are discussed. In [3 , one also 



finds a detailed discussion of the mathematical properties 
of yplV) as well as numerous applications, but this time 
emphasis is on electronic structure theory. 

A new differential operator is not necessarily a useful 
thing, let alone a major achievement. As remarked above, 
differentiating an irreducible spherical tensor F™{r) of the 
type of (jl.ip with respect to the Cartesian components of 
r produces messy expressions. So, if we look at (I8.3P and 
take into account (|8.4p . we may get the impression that 
we replaced something complicated - the left-hand side of 
(|8.3p - by something even more complicated. 

We thus arrive at the paradoxical conclusion that the 
differential operators ^^(V) and 3^™(^) ^'"^ practically 
useful only if it is not necessary to differentiate an irre- 
ducible spherical tensor F^{r) with respect to the Carte- 
sian components of r via the defining explicit expression 
(18. 3p and (18. 4p . Fortunately, this is the case. Differentia- 
tions with respect to Cartesian components can be avoided 
completely. The Laplacian is an irreducible spherical 
tensor of rank zero, and y™(V) is just like the correspond- 
ing regular solid harmonic y"^{r) an irreducible spherical 
tensor of rank £ (a formal proof can be found in [6] p. 
312]). Consequently, we can hope for substantial computa- 
tional and technical benefits if products involving y]^{V) 
and other irreducible spherical tensors arc handled via the 
powerful machinery of angular momentum coupling. 

The first article on the differential operator 3^™(V), 
which I am aware of, is due to Hobson who derived in 
1892 a very consequential theorem on the differentiation 
of functions /: M" C [2311 p. 67]. This theorem is also 
discussed in Hobson's book 1227', pp. 124 - 129], in Martin's 
book 4, Section 3.5], and in 3, Section 3]. 

With the help of Hobson's theorem, it is compara- 
tively easy to obtain explicit expressions for the product of 
3^™(V) and a radially symmetric function (p: M.^ ^ C that 
only depends on r = |r| (see for example |3i Eq. (3.4)]): 



yr(^)v{r) 



1 d 

r dr 



ip{r) 



yrir). 



(8.5) 



Thus, 3^™(V) is a generating differential operator that 
transforms an irreducible spherical tensor of rank zero to 
an irreducible spherical tensor of rank £. 

With the help of Hobson's theorem it is also 
possible to construct explicit expressions for products 
y™^ {\/)F^^ {r). For that purpose, let us assume that the 
irreducible spherical tensor F^^{r) satisfies a relation of 
the kind of (|8.5p . i.e., it can be generated by applying 
3^™^(^) to a suitable scalar function ^^^(r): 



^rw = 3^r(v)ci>,,(r). 



(8.6) 



In view of = yr^yr^-M'^^^ir), we 

need an explicit expression of manageable complexity for 
the product 3^"HV)y^"' (V). This can be accomplished 
easily. Since Gaunt coefficients defined by (jA.6p linearize 
the product of two spherical harmonics, multiplication of 
(|A.7p by yields the linearization formula for the 
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product of the regular solid harmonics (r) and Y^™^ (r) 
(see for example 3, Eq. (3.6)]). In this linearization for- 
mula, we only have to replace the Cartesian components 
of r by those of V and obtain (see for example 3, Eq. 
(3.7)]): 

^max 



^ ^2Aeyp+,n,^^y (8.7) 



The abbreviation Ai, which is either a positive integer or 
zero, is defined by (|A.9|) . 

By combining ()8.5|) . (j8.6p . and (|8.7p we obtain (see for 
example [3l Eq. (3.8)]): 

fmax 

(V) (r) = £ (^mi + m2 |^imi l^ama) 

3;fi+'"^(r). (8.8) 



2A£ 



r dr ' 



X V 



In principle, ()8.8p should suffice for our purposes since 
the scalar function <&£2(r) in (|8.6p can according to (|8.5p 
be obtained from the scalar function fi^ (r) in by re- 
peated integration with respect to r. However, repeated 
integrations can at least potentially lead to nontrivial tech- 
nical problems. Alternative expressions for the product 
yp^^ (V)F™^ (r) are thus deshable. 

By systematically exploiting the tensorial nature of 
the spherical tensor gradient operator, the product 
(V)F^™^ (r) can be expressed as a finite linear combi- 
nation of Gaunt coefficients, radial functions ^fj^i^ir), and 
spherical harmonics [2321 Eq. (4.7)]: 

yr:{^)F^:^{r) 

= ^ ^'^^imi + m2\eimi\t2m2) 



(8.9) 



The functions 'Ji-^g^ {r) in (|8.9p can be obtained by differen- 
tiating the radial part f^^ {r) of the spherical tensor F™^ (r) 
with respect to r = \r\ (see Eqs. (3.29), (4.15) - 

(4.18), and (4.24)] or [21 Eqs. (4.11) - (4.16)]): 



i M i-AiU-a{i)-l/2), 

^1^2 in - 2^ -] ^ r 

q=o ^• 



1 d 

r dr 



(8.10) 



1 d 
r dr 



Al 



r dr 



(8.11) 



1 d 



AI2 



l_d 

r dr 



Ai 



(8.12) 



1 d 
r dr 



A£2 



1 d 

r dr 



A£2 



X r 



-2£-l 



1 d 

r dr 



t2 + l 



Mr) 



.13) 



l_d 

r dr 



«i+^2-3£-l 



l_d 

r dr 



A£ 



A£2 

E 



r dr 



(-A^2).(A4 + l/2). ,, 
; z r 



.14) 



-2s-l 



1 d 

r dr 



£2+1 



(8.15) 



The abbreviations Al, A/i, A/2, and (j{£) are defined by 

(IX9I1 - (KM- 

Other expressions for the product y^^ {V)F^ (r) can 
be found in articles by Bayman |228| . Santos [23lj . Stuart 
], Niukkanen |235| . and Rashid 



There are some radially symmetric functions of con- 
siderable relevance in electronic structure theory that lead 
to remarkably simple expressions if yj^(V) is applied to 
them via (|8.5p . The classic example is the Coulomb po- 
tential 1/r. Hobson I231j showed that the irregular solid 
harmonic Zp{r) — r~^~^Yf^{6,(l)) is generated by apply- 
ing yp{\7) to l/r (further details can be found in Hobson's 
book |2271 pp. 124 - 129]). In modern notation, Hobson's 
result can be expressed as follows (see for example [1211 
Eq. (4.16)]): 



.16) 



Thus, (|8.3p and (|8.16p imply that Cartesian derivatives of 
the Coulomb potential as considered by Guseinov in [321 
Section 3] or by Guseinov and Mamedov [60] can be ex- 
pressed by linear combinations of irregular solid harmonics 
multiplied by integral powers of the Laplacian V^. If 
a nonzero power of the Laplacian acts on an irregular solid 
harmonic, we obtain (see for example [2371 Eq. (29)]): 



.17) 



Here, ^™ is the so-called spherical delta function (see for 
example [2371 Eq. (30)]): 



STir) 



i-lY 



(2£- 1)!! 



3^r(V)<5(r) 



(8.18) 
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If we set in ()8.17|) £ = 0, we obtain tlie well known Poisson 
equation of a unit point charge: 



V^- = -ATTS{r) 



.19) 



Thus, the spherical delta function (5™ can be viewed as 
a generalized solution of the Poisson equation of a unit 
multipole charge. 

In [SlISllMISaiMllIlllllllllSlllSllSSllSTlEI], 

Guseinov considered Cartesian derivatives of the Yukawa 
potential and more generally of Slater- type functions. Also 
in these cases, is is advantageous to rewrite differentiations 
with respect to Cartesian components in tensorial form ac- 
cording to ()8.3|) . 

In this context, it is highly advantageous that Slater- 
type functions with integral principal quantum numbers 
can according to (jS.lSp be expressed as finite linear combi- 
nations of B functions. The reason is that it is remarkably 
easy to apply the spherical tensor gradient operator to B 
functions. For example, application of y^{y) to a scalar 
B function simply produces a nonscalar B function |121|, 
Eq. (4.12)]: 



i-aY 



20) 



It is also fairly easy to apply the spherical tensor gradient 
operator to a nonscalar B function 232, Eq. (6.25)]: 

fmax 



^ /AA 

E(-l)*(t ]BZ't,:-,-tMr). (8.21) 



If we set n = ^ = m = in (|5.15p . we see that the 
Yukawa potential is also a special B function: 



Pk-miM = (47r)i/2/3Bo%(/3,r). (8.22) 



Thus, (|8.20p implies that the application of the spherical 
tensor gradient operator to the Yukawa potential yields the 
so-called modified Helmholtz harmonic |219[ Eq. (6.9)]: 

B!",,,(/3,r) = {AnY'\-l3r'yT{V)Bl,{l3,r). (8.23) 

In the limit of vanishing screening {(3 — > 0), the modi- 
fied Helmholtz harmonic approaches an irregular solid har- 
monic according to |143l Eq. (3.10)] 

ZT{r) = m - l)!!]-i hm [/3^+ii?LV(/3, r)] . (8.24) 

Moreover, the modified Helmholtz harmonic satisfies the 
following exponentially screened variant of the anisotropic 



Poisson equation ((8T7| [3, Eq. (6.18)]: 
[l-r'V^] i?!!'«(/3,r) 



(-1/ 
47r 



An 



(2£-l)!!<5r(r). 



.25) 



In view of this relationship and also because of [2321 Eq. 
(5.6)] 

[1 - V^] BZ]e{P, r) = i?^i,,(/3, r) , (8.26) 

which shows that the differential operator 1 — of 
the modified Helmholtz equation acts as a ladder opera- 
tor in the case of B functions, it makes sense to define a 
distributional B function as the following derivative of the 
three-dimensional Dirac delta function [232|, Eq. (6.20)]: 

(2^- l)!!47r 



B' 



k-e,e 



fc e N. 



.27) 



It is also easy to apply an integral power of the Lapla- 
cian to a i? function. The binomial expansion of (3~'^'^ V^"^ 
in powers of 1 — l3~'^W'^ in combination with (|8.26p yields 
[232, Eq. (5.7)]: 



Em: 



.28) 



The B function relationships given above show that the 
application of the spherical tensor gradient operator mul- 
tiplied by integral powers of the Laplacian to the Yukawa 
potential leads to remarkably compact expressions. Since 
(|8.28p also holds for distributional B functions of the type 
of (|8.27p . it is almost trivially simple to keep track of delta 
function contributions. 

If only first and second order derivatives with respect to 
the Cartesian components of the Coulomb or the Yukawa 
potential are needed, it is not really necessary to do the 
differentiations via (|8.3p . and it is also not too difficult 
to keep track of delta function contributions. If, however, 
differentiations of (very) high orders have to be done, the 
tensorial version (|8.3p offers substantial advantages. 

For example, the two-range addition theorems dis- 
cussed in [T1I2] and in [3 Section 7] were derived via the fol- 
lowing decomposition of the translation operator in (|1.3p 
in terms of tensorial invariants: 



„2k y2fc 



< 



> 



oo 

2"^+2'=fc!(l/2),+fe+i 



.29) 



Apparently, this expansion was first published by Santos 
[233, Eq. (A. 6)], who emphasized that it should be useful 
for the derivation of addition theorems, although he never 
used it for that purpose. However, the addition theorems 
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derived in [U [2] and in 'W, Section 7] show that the ten- 
sorial decomposition (|8.29p is indeed a practicaUy useful 
mathematical tool, although it involves irreducible spher- 
ical tensors of arbitrary order. 

After this essentially technical digression on Guseinov's 
way of differentiating functions /: R'^ ^ C with respect 
to the Cartesian components of the argument vector, more 
fundamental questions such as the existence and conver- 
gence of derivatives of one-range addition theorems will be 
discussed. 

As emphasized already several times, two-range addi- 
tion theorems are essentially rearranged three-dimensional 
Taylor expansions in the Cartesian components of the shift 
vector, which converge pointwise and in suitable open sets 
even uniformly. Accordingly, two-range addition theorems 
represent ordinary functions, which are locally smooth and 
can be differentiated. Thus, it is possible to generate more 
complicated two-range addition theorems by differentiat- 
ing simpler ones (see for example [1161 1219] ). 

In contrast, Guseinov's one-range addition theorems 
discussed in Sections H] and [7] - if they exist at all - either 
converge in the mean with respect to the norm of an ap- 
propriate Hilbert space or they converge weakly, i.e., they 
only make sense in suitable functionals restricted to sub- 
sets of certain Hilbert spaces. As is well known, smooth- 
ness is neither necessary for strong or weak convergence 
nor implied by it. Thus, it is not clear whether Guseinov's 
one-range addition theorems represent sufficiently smooth 
functions that can be differentiated, and if yes, in which 
sense they can be differentiated. 

These problems are not uncommon in mathematics. 
Often, it is comparatively easy to prove convergence in the 
mean, but very difficult or even impossible to prove point- 
wise convergence or even uniform convergence directly. 
Scenarios like this one obviously raise the question, under 
which conditions uniform convergence is implied by con- 
vergence in the mean. A condensed review of these issues 
and a discussion of some elementary results can be found 
in a recent article by Ford and Pennline [238j . Related 
topics, albeit with an emphasis on distribution theory, are 
discussed in the book by Strichartz ^239, Chapter 8]. 

Guseinov's derivatives of strongly convergent one-range 
addition theorems for Slater-type functions can be made 
mathematically meaningful if it can be shown that these 
addition theorems do not only converge in the mean, but 
also pointwise and uniformly. 

Probably, there is some hope in the case of the one- 
range addition theorems (j6.ip . which are expansions in 
terms of Guseinov's functions k'^l^£{j,t') and which con- 
verge with respect to the norm (j4.19p of the weighted 
Hilbert space L^fc(R^), if the principal quantum number 
TV of the Slater-type function hiP^'"' ^ ''"') ^ posi- 
tive integer e N. We would have to check whether the 
additional conditions discussed in Szego's book [1151 The- 
orem 9.1.5 on p. 246] are satisfied in this case. Such an 
approach may also be successful in the case of the deriva- 



tives of the one-range addition theorems for Guseinov's 
functions r) considered in 

I am, however, very skeptical about all addition theo- 
rems for Slater-type functions ±r') with nonin- 
tegral principal quantum number A^ G M \ N. 

The situation is more complicated if weakly convergent 
one-range addition theorems of the type of (|7.2p are to be 
differentiated. Weakly convergent addition theorems are 
essentially expansions of generalized functions in the sense 
of Schwartz |129j . whose derivatives also have to be in- 
terpreted in this sense. If / is a generalized function and 
^y; : M'^ ^ C is a suitable test function, derivatives of / are 
defined by the following functional [129, Chapter II]: 

/[(s)"(l)'(i^)'H*'^'^'^ 

(8.30) 

If ^ e 5(R3) or i/i G P(R3), V is infinitely differentiable 
and it is possible to define weak derivatives of arbitrary 
order of a generalized function / via (|8.30p . 

There is a detailed mathematical literature on deriva- 
tives of generalized functions and their interpretation and 
application. For example, there is a monograph by Ziemer 
on weakly differentiable functions ^240 . Unfortunately, 
the mathematical literature is not particularly helpful if 
we want to apply derivatives of one-range addition the- 
orems in multicenter integrals and if we insist on using 
exponentially decaying functions. Mathematicians almost 
exclusively use infinitely differentiable test functions be- 
longing to either 5(M'^) or I?(R'^). However, Kato's work 
on cusps [1?5^ shows that atomic and molecular wave func- 
tions are not infinitely differentiable. This is also true if we 
approximate atomic and molecular wave functions in vari- 
ational calculations by exponentially decaying functions. 

If we want to use one-range addition theorems, which 
are derivatives of generalized functions, in multicenter in- 
tegrals, we are confronted with the annoying problem that 
the remainder of the integrand, which assumes the role of 
the test function, not only has to decay sufficiently rapidly 
for large arguments, but also has to possess continuous 
derivatives of a sufficiently high order. Obviously, this 
complicates considerably our attempts to prove that the 
resulting expansions converge to the correct results. Math- 
ematicians know why they prefer infinitely differentiable 
test functions belonging to either iS(R'^) or X'(R^). 

To illuminate these problems, let us apply the spherical 
tensor gradient operator yf^^iy) — yfjiSr) to the one- 
range addition theorems (|7.2p for the Coulomb potential. 
Then, we formally obtain with the help of (|8.16p one-range 
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addition theorems for the irregular sohd harmonic: 



{2L-1)\ 



r — r' 



— ll!! ^ ^ n'i'm' \P ) 

n' £'m' 

X [3^f(V) fe*™ (/3,r)] ,vl;;5>(/3,r'). 



.31) 



An essentially identical expression - the factor (—1)'^ 
would be missing - can be derived by differentiating 
l/\r — r'\ instead with respect to r'. 

If we use the addition theorems (|8.3H) in a multicen- 
ter integral, we essentially have two options: Either, we 
can try to construct an explicit expression for the prod- 
uct y^(y) k'i";^i{P,r), or we could use (fg:^ and apply 

{V) to the r-dependent part of the remainder of the 
integrand, which assumes the role of the test function. 

If we apply J^^^(V) to fe^'™^(/3,'r) and use one of 
the numerous expressions for Je^e^ir) defined by jHUI), it 
should be possible to construct a closed form expression 
for this product. However, the resulting functions of r will 
no longer be orthogonal and they may also have a (much) 
more complicated structure. 

Alternatively, we can use (|8.30p and shift 3^|^(V) to 
the remainder to the r-dependent part of the integrand. 
If this happens to be a nonclassical two-center density, we 
have to apply the Leibniz theorem of the spherical ten- 
sor gradient operator (see jj] Eq. (4.22)] and references 
therein). Messy expressions are then likely, possibly in- 
volving singular or distributional contributions. In either 
case, nontrivial technical difficulties are to be expected. 
One must not forget that we still have to prove that the 
resulting series expansion for the multicenter integral con- 
verges to the correct result. 

If we should indeed need weakly convergent one-range 
addition theorems of the type of ()7.2p for the irregular 
solid harmonic Zff {r — r'), then it is probably easier to 
construct it from the scratch by expanding Z^^ in terms of 
Guseinov's functions fc^'™£(/5, »")• I strongly suspect that 
in this case it would be much easier to prove convergence 
of the resulting expansion to the correct result. 

One the basis of our current level of knowledge (or 
rather the lack of it), the use of derivatives of weakly con- 
vergent addition theorems in multicenter integrals would 
be purely experimental. 

9 Summary and Conclusions 

The efficient and reliable evaluation of multicenter inte- 
grals is among the oldest mathematical and computational 
problems of molecular electronic structure theory (a re- 
view of the older literature can be found in an article by 
Dalgarno [241j ). In spite of the efforts of numerous re- 
searchers including Guseinov and coworkers, no completely 



satisfactory solution has been found yet. The situation is 
particularly unsatisfactory in the case of the notoriously 
difficult integrals of exponentially decaying functions, but 
even in the case of the much simpler integrals of Gaus- 
sians, there is still active research going on (see for example 
[ini EUl Em ETH eh eh Eig and references therein) . 

Different centers in the integrand make the evaluation 
of an integral difficult. They prevent the straightforward 
separation of a multicenter into products of simpler inte- 
grals. As discussed in Section[l] addition theorems, which 
are expansions of a given function /(r±r') with r, r' G R'^ 
in products of other functions depending on either r or 
r', are principal tools that can accomplish a separation of 
variables, albeit at the cost of infinite series expansions. 

In electronic structure calculations, predominantly 
those addition theorems have played a major role that 
depend on r and r' only indirectly via r< and r> and 
thus possess a two-range form. The prototype of such 
an addition theorem is the Laplace expansion (|1.2p of the 
Coulomb potential. Unfortunately, the use of two-range 
addition theorems in multicenter integrals can easily lead 
to nontrivial technical problems. This explains why many 
authors have tried to construct one-range addition theo- 
rems that can be applied more easily. 

Two-range addition theorems can be constructed by 
applying the translation operator e''<'^> in its tensorially 
invariant form ()8.29p to /(t">), which is by assumption 
analytic at r>. Accordingly, two-range addition theorems 
are essentially rearranged three-dimensional Taylor expan- 
sions in the Cartesian components of the shift vector r< 
that converge pointwise and in suitable open sets even uni- 
formly. If / is not analytic at the origin, its addition the- 
orem converges only if it has a two-range form. 

Classical and complex analysis is mainly concerned 
with power series that converge pointwise. However, point- 
wise convergence is a very demanding requirement, and its 
scope is too limited to cover all cases of interest. If we 
want to avoid the troublesome two-range form of addition 
theorems, we have to replace pointwise convergence by a 
weaker form of convergence. 

One-range addition theorems can be constructed by 
exploiting Hilbert space theory, whose basic facts are re- 
viewed in SectionEl In this approach, a function /(r ±r') 
belonging to a suitable Hilbert space is according to p.ip 
expanded in terms of a function set {(p™i{r)}n,e,m that is 
complete and orthonormal in this Hilbert space. In gen- 
eral, such an expansion converges only in the mean, i.e., 
with respect to the norm of the underlying Hilbert space, 
but not pointwise. 

The natural Hilbert space for electronic structure cal- 
culations based on effective one-particle wave functions is 
the Hilbert space L^(M'^) of square integrable functions de- 
fined by (|2.8p . As discussed in Section [31 it is also possible 
to construct one-range addition theorems that converge 
with respect to the norm of a suitable weighted Hilbert 
space Ll{R^) defined by ([2ll|l . Here, w{r) > is a suit- 
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able weight function. There is, however, a principle prob- 
lem: In general, we neither have L^(R'^) C i^(R'^) nor 
L^(K3) C L2(M3). Accordingly, it is by no means easy 
to find a nontrivial weight function 'w{r) ^ 1 that leads 
to substantial computational benefits and simultaneously 
avoids major drawbacks. 

If we demand that the complete and orthonormal func- 
tions {'p™i{r)}n,e,m are also irreducible spherical tensors 
of the type of (|l.ip . we more or less automatically arrive 
at functions based on the generalized Laguerre polynomi- 
als. Section m discusses Lambda functions A™^(/3, r) de- 
fined by (|49|l . Sturmians 4'™^(/3,r) defined by ([?1T|1 . and 
Guseinov's functions fc^'™f(/?,'r) with k = —1,0,1,2,... 
defined by (gH]). 

Guseinov's functions satisfy the orthonormality rela- 
tionship ()4.17p involving the weight function w{r) — . 
Accordingly, they are complete and orthonormal in the 
weighted Hilbert space ^^^(K^) with k = -1,0,1,2,... 
defined by (j4.20p . Guseinov's functions provide some uni- 
fication since they contain for k — Q and fc = — 1, respec- 
tively. Lambda functions and Sturmians as special cases. 
However, this unification does not apply to approximation 
processes. For fc / 0, we neither have ^^^(R'^) C L-^(R^) 
nor i2(R3) ^ Ll^{m.^). Thus, / £ L'^{^^) does not imply 
that an expansion of / in terms of Guseinov's functions 
with fc ^ necessarily converges in L'^^{R^). 

One-range addition theorems for exponentially decay- 
ing functions had been considered by several researchers 
before Guseinov. The so far most compact results were 
obtained by Filter and Steinborn [T^ who derived the 
symmetrical one-range addition theorems (|5.4p and (|5.5p 
for Lambda and B functions, respectively, by expanding 
them in terms of Lambda functions, which are complete 
and orthonormal in L^(R'^). In Section [5l it is shown that 
the approach of Filter and Steinborn, which was based on 
the convolution theorem (|5.6p of B functions, works also 
in case of Slater- type functions. The approach of Filter 
and Steinborn can be generalized to expansions in terms 
of Guseinov's functions. However, these addition theorems 
have a more complicated structure than the corresponding 
expansions in terms of Lambda functions. 

Section [S] describes Guseinov's mathematically dubi- 
ous attempts of deriving one-range addition theorems for 
Slater- type functions. Guseinov first derived the expan- 
sions ()6.ip of a Slater-type function Xn hiP'''' ^ with 
integral or nonintegral principal quantum number N in 
terms of his functions fc^™^(7,f") with 7 > 0. Guseinov's 
results are mathematically correct, although not necessar- 
ily optimal. In a later step, Guseinov replaced his or- 
thonormal functions by nonorthogonal Slater-type func- 
tions with integral principal quantum numbers according 
to (j6.4p and rearranged the order of summations in his 
addition theorems. In this way, Guseinov obtained the ex- 
pansions (|6.8p of hiPi''' ^ ™^ terms of Slater- type 
functions X™f(7i with integral principal quantum num- 
bers n e N. 



As reviewed in Section [21 a function / belonging to a 
given Hilbert space can be expanded in terms of complete 
and orthogonal functions. Such an expansion converges 
with respect to the norm of the Hilbert space, and the co- 
efficients of this expansion satisfy Parseval's equality (|2.5p . 
Thus, the expansion coefficients are bounded in magnitude 
and they also vanish for large indices. 

If the orthogonal functions are replaced by nonorthog- 
onal functions, it is possible to construct finite approxima- 
tions of the type of (|2.2p . Unfortunately, it is not guaran- 
teed that an infinite expansion of the type of (|2.4p exists, 
and if it does, it can happen that the coefficients of this 
expansion are unbounded and do not vanish for large in- 
dices. 

Accordingly, Guseinov's manipulations, which pro- 
duced the rearranged addition theorem (|6.8p from (j6.ip . 
are dangerous. Thus, he should have proved that his ma- 
nipulations are legitimate and lead to a meaningful result. 

Addition theorems of exponentially decaying functions 
are fairly complicated mathematical objects. Therefore, 
it is very difficult or practically even impossible to prove 
explicitly that a given addition theorem converges or di- 
verges. This applies also to Guseinov's rearranged addi- 
tion theorem (|6.8p for Slater-type functions Xn hif^^ r±r'). 
However, as shown in Section[Sl the one-center limit r' = 
of Guseinov's addition theorem (|6.8p does not exist if the 
principal quantum N is not a positive integer. Thus, for 
iV S R \ N the rearranged addition theorem ()6.8p does not 
exist for the whole argument set R"^ x R''. The remaining 
question, whether the rearranged addition theorem (|6.8p 
exists for iV e N, cannot be decided in this way and is still 
open. 

The nonexisting rearranged addition theorems (|6.8p 
with N — L — M = were the starting point for 
Guseinov's construction [HH] of one-range addition theo- 
rems for the Coulomb potential, but this is not the only 
problem. The Coulomb potential does not belong to any 
of the Hilbert spaces implicitly used by Guseinov, since 
they all involve an integration over the whole R'^. Ac- 
cordingly, expansions of l/|r — r'| in terms of Guseinov's 
functions r), which formally yield the symmet- 

rical one-range addition theorems (|7.2p . diverge for any 
fc = —1, 0, 1, 2, . . . in the mean with respect to the norm 
of the weighted Hilbert space L^fc(R'^). 

This observation seems to imply that all attempts of 
constructing a symmetrical one-range addition theorem for 
l/|r — r'l by expanding it in terms of functions, that are 
complete and orthonormal with respect to an inner prod- 
uct involving an integration over the whole M"^, are futile. 
Although seemingly obvious, this conclusion is premature 
and the situation is better than it may look. 

It is possible to use divergent one-range addition theo- 
rems like (j7.2p in multicenter integrals in a mathematically 
rigorous way. However, convergence in the mean, which in 
the mathematical literature is frequently called strong con- 
vergence, is too demanding for that. We have to replace 
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it by an even weaker form of convergence. Instead, we re- 
quire that a possibly divergent one-range addition theorem 
produces meaningful results when used in a multicenter in- 
tegral. Obviously, this approach is inspired by the theory 
of generalized functions in the sense of Schwartz |129j . 

Weak convergence of addition theorems is all we really 
need in variational electronic structure calculations. Point- 
wise convergence or convergence in the mean - although 
undeniably convenient - are not really necessary. Thus, 
the use of weakly convergent addition theorems offers new 
perspectives and simplifies or even solves certain technical 
problems. However, it would be overly optimistic to ex- 
pect that we could get all these benefits without having to 
pay a price. 

The key problem, which makes the application of 
weakly convergent addition theorems difficult, is the for- 
mulation of sufficiently simple criteria which guarantee the 
convergence of the resulting expansions for multicenter in- 
tegrals. In that respect, a lot of work remains to be done. 

Regularity conditions - or rather their absence - also 
play a major role in Section[Sl In several articles, Guseinov 
had constructed more complicated one-range addition the- 
orems by differentiating simpler one-range addition theo- 
rems with respect to the Cartesian components of their ar- 
gument vectors. Guseinov's one-range addition theorems 
- if they exist at all - either converge in the mean with re- 
spect to the norm of an appropriate weighted Hilbert space 
or they converge weakly, i.e., they only make sense as gen- 
eralized functions in suitable functionals. It is thus not at 
all clear whether Guseinov's addition represent functions 
in the ordinary sense that are locally smooth and can be 
differentiated. We need to know much more before we can 
safely apply derivatives of one-range addition theorems in 
multicenter integrals. 

Section [S] also discusses a more convenient alternative 
to the troublesome differentiation of irreducible spherical 
tensors of the type of (|l.ip with respect to the Cartesian 
components of their argument vectors. The alternative 
proposed here is based on the reformulation (|8.3p of a dif- 
ferential operator as a finite sum of integral powers of the 
Laplacian multiplied by irreducible spherical tensor gradi- 
ent operators 3^™(V) defined by (|8.4p . This differential 
operator plays a decisive role in the derivation of two- 
range addition theorems as rearranged Taylor expansions 

At first sight, the use of (|8.3p looks like a bad idea. 
However, y™(V) is an irreducible spherical tensor of rank 
i. Consequently, products 3^™^ (V)F™^ (r) can be simpli- 
fied considerably with the help of the powerful machinery 
of angular momentum coupling, yielding the finite linear 
combination (|8.9p of Gaunt coefficients, radial functions, 
and spherical harmonics. Thus, no differentiations with 
respect to the Cartesian components of r are needed. We 
only have to differentiate with respect to the radial vari- 
able r = \r\. 

It is the intention of this article to demonstrate that 



the analytical tools of classical and complex analysis do 
not suffice for a successful treatment of one-range addition 
theorems. More modern and also more sophisticated con- 
cepts such as Hilbert spaces and generalized functions are 
indispensable. 

In particular, we need different concepts of conver- 
gence. If a function / belongs to a suitable Hilbert space, 
one-range addition theorems can be constructed that con- 
verge in the mean with respect to the norm of that Hilbert 
space, even if / is not analytic everywhere in M"^ . If, how- 
ever, / does mot belong to that Hilbert space, one-range 
addition theorems can only converge weakly in the sense 
of generalized functions. 

Analogous considerations apply also with respect to 
numerical techniques. In the context of multicenter in- 
tegrals, addition theorems are essentially mathematical 
recipes that generate infinite series expansions. We have 
no a priori reason to assume that these series expansions 
converge rapidly. Thus, our assessment of the usefulness 
of addition theorems depends crucially on our ability of 
evaluating slowly convergent or even divergent series ef- 
fectively and reliably. 

The conventional process of adding up the terms of 
an infinite series successively until convergence is finally 
achieved is in far too many cases hopelessly inefficient. It 
is thus an obvious idea to use convergence acceleration 
and summation techniques. Of course, we do not know 
whether the currently known convergence acceleration and 
summation techniques, which were mentioned shortly in 
Section [71 are powerful enough to produce satisfactory re- 
sults or whether it will be necessary to construct new and 
hopefully more powerful transformations. This has to be 
investigated. 

Let me emphasize: If we want to evaluate even difficult 
multicenter integrals effectively and reliably via series ex- 
pansion, then there is no alternative to the systematic and 
intelligent use of convergence acceleration and summation 
techniques. 

It is dangerous to rely too much on completeness. The 
completeness of a function set in a Hilbert space only im- 
plies that any / belonging to this Hilbert space can be 
expanded in terms of these functions and that this expan- 
sion converges in the mean. Unfortunately, convergence is 
a very weak statement. In particular, the convergence of 
an infinite series does not imply that this series is numeri- 
cally useful. A simple and yet striking example is provided 
by the Dirichlet series C,{s) = X]^^o(" + ^)^^ ^"-"^ Wie- 
mann zeta function. This series converges if Re(s) > 1, 
but is notorious for extremely slow convergence if Re(s) is 
only slightly larger than 1. For example, in }207[ p. 194] 
it was estimated that in the order of 10®°" terms of the 
Dirichlet series would be needed to compute C(l.Ol) with 
an accuracy of 6 decimal digits. 
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A Spherical Harmonics and Gaunt 
Coefficients 

If we choose the phase convention of Condon and Short- 
ley |106l Chapter III. 4], the spherical harmonic can be 
expressed as follows [SJ p. 69]: 



_ jm+|m| 



(2^+l)(^- |m|)! 
47r(^+ |m|)! 

X pl"l(cos6l)e'™'^. 



1/2 



(A.l) 



Here, pj'"'(cos6') is an associated Legendre polynomial 
|1071 p. 155]. Alternative phase conventions for the spher- 
ical harmonics are discussed in [243', pp. 17 - 22]. 



The spherical harmonics Yl 



are often called sur- 



face harmonics because the angles Q and characterize 
a point r/r on the surface of the three-dimensional unit 
sphere. In the literature, it is common to introduce the 
so-called regular and irregular solid harmonics 



ZT{r) = r-'-^Yne^). 



(A.2) 
(A.3) 



The regular solid harmonic y^^ir) is a homogeneous har- 
monic polynomial of degree I in the Cartesian components 
of r = [x, y, z) p. 71]: 



ynr) 



E 

fe>0 



2i+l 
An 

{-X- 



{t + m)\(£ ~ m)\ 



1/2 



\ni\^ 

) 



2™+2'=(m + k)\kl{e - m - 2k)\ 



(A.4) 



Moreover, 3^™(r) is for all r e M'^ a solution of the homo- 
geneous three-dimensional Laplace equation 



dx'^ 



d 



2 1 



fir) 



0, 



(A.5) 





whereas the irregular solid harmonics are generalized so- 
lutions of the Poisson equation (|8.17p of a unit multipole 
charge, yielding the spherical delta function (|8.18p . 

The so-called Gaunt coefficient | 244| is the integral of 
the product of three spherical harmonics over the surface 
of the unit sphere in R'^: 



(4TO3K2m2|^lTOl) 



(A.6) 



It follows from the orthonormality of the spherical har- 
monics that Gaunt coefficients linearize the product of two 
spherical harmonics: 



(2) 



(^mi +m2|£imi|^2m2)r;"^+"'^(r!). (A.7) 



The symbol X)''^'' indicates that the summation proceeds 
in steps of two. The summation limits in (jA.7[) are given 
by [5111 Eq. (3.1)] 



A. 




IS even , 
1 is odd , 



(A.Sa) 
(A.8b) 
(A.Sc) 



A compact review of the properties of Gaunt coefficients 
and additional references can be found in [31 Appendix C] . 
In this article, the following abbreviations are used: 



Ml = (£-4+4)/2, 

M2 = (^ + ^l-^2)/2, 

a{l) = (4+^2+^)/2. 



(A.9) 
(A.IO) 
(A.ll) 
(A.12) 



If the three orbital angular momentum quantum numbers 
li, £2, and £ satisfy the summation limits (|A.8[) . then these 
quantities are either positive integers or zero. 
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